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Laughlin's wave functions, Coulomb gases and 
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, y^. , In the context of the fractional quantum Hah effect, we investigate Laughhn's cele- 

brated ansatz for the ground state wave function at fractional filling of the lowest Landau 
level. Interpreting its normalization in terms of a one component plasma, we find the effect 
of an additional quadrupolar field on the free energy, and derive estimates for the ther- 
modynamically equivalent spherical plasma. In a second part, we present various methods 
for expanding the wave function in terms of Slater determinants, and obtain sum rules for 
the coefficients. We also address the apparently simpler question of counting the number 
of such Slater states using the theory of integral polytopes. 
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1. Introduction 

These notes originate from an attempt to understand the normalization - and other 
properties - of the many body fermionic wave functions suggested by Laughlin as candi- 
dates for the ground states of the fractional quantum Hall effect. Similar quantities appear 
in the context of the two-dimensional classical one-component plasma (sometimes called 
jellium) and in statistical problems related to matrices or random polynomials. 

Our work is divided into three parts. In the first one, after recalling the simplest cases 
of quantum wave functions at odd fractional filling of the lowest Landau level proposed 
by Laughlin, we discuss their expansion in terms of Slater determinants which would 
allow to obtain- among other things - their normalization (section 2). We then record the 
standard interpretation of the normalization problem as computing the partition function 
of a two dimensional - neutralized - one component plasma. In a variant we show that an 
added quadrupolar compensating field has no effect on the thermodynamic properties (at 
complete filling) except to modify the shape of the " quantum liquid drop" . It is possible to 
substitute to the planar geometry a spherical one with higher symmetry (section 3). This 
permits one to obtain bounds for the residual free energy - in the Coulomb interpretation- 
using Holder's and Hadamard's inequalities. These bounds can be conjecturally improved 
by using the behaviour at zero temperature of some special symmetric configurations for 
finitely many particles (up to 42). We devote appendix C to an interesting instability of 
the cuboctahedral configuration. 

Part two is centered around one of the most fundamental objects of algebra: the 
discriminant of a polynomial 

D^l[ix,-x,), (1.1) 

where i runs over A^ distinct values (conveniently chosen as the integers 0, 1, 2, ..., A^ — 1 ) 
and Xi are indeterminates which may be thought of as the roots of the monic polynomial 

N 
i k=0 

The discriminant D as well as its powers are symmetric polynomial functions of the in- 
determinates Xi and as such admit a unique polynomial representation in terms of the 
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elementary symmetric functions a^, a well known fact since the days of Newton. The 
subsequent representation for D is unfortunately not very useful for our purposes. On the 
other hand the ring of symmetric polynomial functions with rational coefficients can be 
viewed as a graded vector space over Q. Each homogeneous subspace admits a basis in 
terms of Jacobi-Schur functions chY{xi), indexed by partitions or Young tableaux. Thus 
with s a positive integer 

(-lym^-^y^D^ix^) = Yl 9^y^ ^^^(^^) • (1-3) 

|y|=sAr(iV-l) 

We were unable to find in the immense literature on symmetric functions and repre- 
sentations of linear groups a discussion of such an expansion. In section 4 we give various 
expressions for the coefficients gy (eqs. (4. 16), (4. 20), (4. 22) and (4.57)) none of them being 
very effective. We also derive a number of their properties and present some tabulations 
in appendix D. As already admitted, even for s = 1, we could not find a useful expression 
for the general integral coefficients gy simple enough for our purposes. This is left as a 
challenge for a clever reader. However in section 5 we exhibit a remarkable sum rule 

\-|>)|2- [i^s + l)N]l 

derived from a formula conjectured by F. Dyson, and proved by Wilson, Gunson, and 
Good. 

The gy coefficients also satisfy an overdetermined linear system (section 6) which 
turns out in practice to be the fastest means of computation for small enough number of 
particles. Finally we consider ^-specializations in section 7. 

We devote the third part to yet another problem, at first sight elementary, which is 
to count the number of terms in the expansion (1.3). It turns out that this is not trivial 
at all. We will elaborate in section 8 and appendix E on our findings in this direction. 

As this work was in progress we received an article by G.W. Dunne |I[] which partly 
overlaps with ours. 

2. From the Laughlin wave function to the classical Coulomb gas 

In the fractional quantum Hall effect one observes plateaux in the Hall conductance 
indexed by filling fractions of Landau levels, the infinitely degenerate non-relativistic en- 
ergy levels for a charged particle interacting with a magnetic field (see for instance for a 



review) . Assuming the particles constrained in a transverse plane, and choosing a suitable 
scale for the field, an orthogonal complete set of states in the lowest Landau level is given 
by 

(2.1) 



ipi{x) = —=x e 
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where a; is a complex variable and / = 0, 1, ..., can be interpreted as the angular momentum 
eigenvalue. The normalization of these wave functions is 
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For a system of N non-interacting particles with Fermi statistics in the lowest Landau level 
we take Slater determinants as an orthogonal basis of wave functions 
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(2.3) 



We may assume ordered indices < /q < /i < ... < /tv-i, since the multiparticle wave 
function is antisymmetric in the permutations of the Vs. We have then 
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{K,...M-Mi'o,-,i'^J = Nll[hlS^ 
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(2.4) 
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The state with lowest total angular momentum, corresponds to the choice /« = i, in which 
case the Slater determinant reduces up to a factor to the Vandermonde determinant 
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(2.5) 



We use the notation D{xo, ..., xn-i) for the discriminant, equal up to a sign to the square 
of the Vandermonde determinant 



D{xo,...,xn- 



N(N-l) 



Y\_ixi - Xj) = (-1) '2 ' A{xo,...,xn-i] 



(2.6) 



and identify the normalization of the "densest state" (filling fraction 1) as a partition 
function through 

N-1 ,2. N 



Zn{1) = / n (^e---)A(x)A(x) = n j! 

■^ i=0 j=l 



(2.7) 



Electrons subject to an intense magnetic field are nevertheless still interacting via 
Coulomb forces and also with the substrate as well as impurities. To this day it is still a 
difficult problem to understand the properties of such a system even at very low tempera- 
ture. Nevertheless the remarkable stability of such quantum fluids exhibited in particular 
by the existence of the Hall plateaux led Laughlin to assume simple trial wave functions at 
fractional fillings. In the simplest cases of filling fractions of the form v = 2jxi 5 s = 0, 1, ... 
(where the odd denominators are required by Fermi statistics), his guess takes the form 



Xo,...,XAr-i 



's V^U? •••? 



1 V — > N — 1 X ■ X ■ 

-p^A(a:o,...,XAr_i)2-+ie-^.-o ^, (2.8) 



such that as Xi -^ Xj the wave functions have zeroes of order 2s + 1. It is then natural to 
ask the following questions: 
(a) compute the normalization 
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Zj,{2s + 1) = (</>,!</>,) = / n (^e--^-0 (a(x)A(x))''^' , (2.9) 

with Z]^{1) = ((/)i|(/)i) given above, 
(b) Find the expansion of Laughlin's wave functions (ps on the complete basis of (free 
particle) Slater determinants 

^s{xq,...,xn-i)= ^ g\^l^^^i^_^^i^^,„^i^_^{xQ,...,XN-i), 

0<Io<...<In-i<{2s+1){N-1) 

(2.10) 
or equivalently 

A(xo,...,XA._i)2^+i= Yl alt.M-. l^''-^'~-M, (2.11) 

0<io<---<^iv-i<(2s+l)(Ar-l) 

with the short hand notation |x^°...x^^-i| for the determinant obtained by replacing 
X successively by xq.xi, ... in the first, second ... line. 
This in turn amounts to 

„ sAf(Af-l) 

A{xo,...,xn-i) ={-l) 2 D{xo,...,XN-iy 

Yl 9if,...,i^_,chi,,„„i^_^ixo,...,XN-i) ^^'-^^^ 

0<Zo<---<'iv-i<(2s+l)(Ar-l) 



The notation c/i... stands for the polynomial characters of the general linear group GL^, 
expressed here for a diagonal matrix X = diag{xo, .., xn-i) (we assume for the time being 
that none of the differences among x's vanishes) as 

chi,,,„^i^_,{xQ,...,XN-i) = 1^^ "^"jv-i| • (2-13) 

Prom L'Hospital's rule this remains meaningful when any subsets of x's coincide. We 
call them general Schur functions - although they were introduced by Jacobi - and shall 
express them in terms of the (normalized) traces 

N-l 

Parenthetically we will also have use for to = N. Defining the sequence 

< /o < /i < ■ ■ ■ < fN-i < 2s{N - 1) (2.15) 

through 

h = n + i i = 0,l,-,A^-l, (2.16) 

one can record the increasing sequence li in a Young tableau with /o boxes in the A^-th line, 
/i boxes in the A^ — 1-th line,..., Jn-i boxes in the first line. We will sometimes switch from 
the notation {li} to y = {fi}. We observe that the expression ( p.l3|) extends naturally 
as an antisymmetric quantity under permutations of /q, ..,/Ar-i and correspondingly for 

The expansion ( |2.10|) allows one to perform the integral in the partition function 

N-l 

Z^{2s + l)=m Yl \9lt.,i.-f n ^'' • (2.17) 

0<Zo<---<;iV-i<(2s + l)(Af-l) i=0 

This formula is useful only insofar we can obtain some mastery of the coefficients gy ■ 



The integral in ( p.9|) continues to make sense for arbitrary (positive) exponent p instead 
of the odd integer 2s + 1, in which case 
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Zn{p) = / n (^e---0 (a(x)A(x))' 

CL Xi — y XiXi-\-2p} ln\ 



'-° (2.18) 

Ig /— /i = "' ' ' '' ^-^0<i<j<N 

i=0 



This suggests two alternative descriptions of our problem. In the first, consider the set 
of monic polynomials of degree N (with complex coefficients) identified with ((D )^y™™ 
and define a probability distribution on this space by considering the roots as independent 
Gaussian variables 

di,N{xo....XN-i) = TT(^^e-^^^0 (2.19) 

i=0 

This is indeed symmetric in any permutation of the x's. The first form of ( |2.18| ) at least 
for p a positive integer involves the computation of the moments of the distribution of the 
module of the discriminant \D\, thus for any positive p 
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The second and more useful interpretation of the integral ( |2.18| ) is as a canonical 
partition function of a one component plasma (logarithmic repulsive two body potential), a 
classical two-dimensional Coulomb fluid, with neutralizing background. In this framework 
one can avail oneself of a considerable body of knowledge in the thermodynamic limit 
A^ — i> OO P] , and in particular of some evidence for a flrst order "fluid-solid" transition in 
the parameter /3 = 2p of the order of 

Arans. = 2ptrans. ^ 142, (2.22) 

according to Choquard and Clerouin [Q]. It would be of great interest to be able to 
understand from flrst principles the mechanism of this transition to a classical Wigner 
solid. 

Following Jancovici and Alastuey, the plasma interpretation justifles a large A^ behav- 
ior for lnZj^{p) obtained as follows. Consider A^ classical particles of charge e in a disk of 
radius R interacting via a repulsive two-dimensional Coulomb potential 

e^ 1 



^ /I U^i ?' 



choosing an arbitrary length scale inside the logarithm (which should anyhow disappear 
from the final result). The contribution of the charges to the total energy is 



e2 



^ch. = -— /n|A(xo, ..,XN-i)\. (2.23) 

47r 



We assume a uniform neutralizing background with density :;;r^- Its self-interaction is 

^back. = f i^ ) ^ / (fxi (fx2 ^lnj^-^—-r (2-24) 

Set 



7ri?V 2 7|^,|,|^2|<i^ 27r |xi-X2| 



(f)(x) = — / d^y In- = H In— \x\ < R. (2.25) 

"^T^ J\y\<R \x-y\ 4 2 R 



This expression follows from the fact that —A(f){x) = 1 for |a:| < -R, while for |a:| = -R, (^ is 
equal to the potential created by the total charge nR^ located at the center, i.e. -^^^In^. 

After a short calculation 

N^e'^ 111 
^back = ^-;^(- + -In — ). (2.26) 

4nR^^8 2 R' ^ ' 

The last contribution to the energy (ignoring kinetic energy which factorizes from the 
partition function) is the attractive interaction between charges and background 

i i 

Thus the total energy divided by /c^T (T the temperature) is 

E 
-, — 77; = {Ech. + -Eback. + Eint.)/kBT 
kbT 

"^ A 3 

(^ Xl^^^ - /n|A(a;o, ..., xn-i)\^ + NHnR - -A^). 



(2.28) 



AnkeT^R^ 
We set 



"" (l=p^ r2 =pN (2.29) 



A-kUbT 2 
in such a way that, 

E Y^ _ ^ , ,., ,, pA2 3pA2 



ksT 



^XiXi -2p ln\A{xo,...,XN-i)\ + -—ln{pN) -— . (2.30) 



The excess free energy -Fexc. of this system, as compared to the perfect gas (i.e. deriving 
only from the kinetic energy), is given by 

^- Inf n^e-^ 



ksT J\xA<R ^^^ 7ri?2 



^=0 (2.31) 

^ -ln{pN) + -pN"^ - Nln{pN) + lnZN{p). 



2 ^^ ' 4^ 

According to thermodynamics, in the hmit of large A^ (i.e. large R) we have extended 
the integrals to the full plane up to errors of order at most yN in -Fexc. (i-e. arising from 
boundary terms). Since we expect an extensive excess free energy, we obtain the estimate 

InZ^ip) = ln{pN) - -pN^ + Nln{pN) + 0{N). (2.32) 

It is the term of order A, namely —-jf^^ depending on p, which is shown to undergo a 
first order phase transition, namely the limit 

1 vN"^ 3 

lim —{lnZN{p)-- ln{pN) + -pN^ - NlnipN)) (2.33) 

N^cxD N 2 4 

should exist and have a discontinuous derivative for a value of p close to 71. In principle 



eqn. (|2.32|) answers question (a) up to the unknown term of order A. But its status is at 
best heuristic. It is therefore reassuring that at least for p = 1 (s = 0) where we have an 
exact result, we can justify (|2.32| ). Note that in the thermodynamic reasoning the value 
of p did not play a crucial role. We devote appendix A to the asymptotic evaluation of 
/nZjv(l) using Euler-Mac Laurin's formula. It reads 

N 

lnZN{l) = lnl[j\ 

^ ^luN - '-N^ + A/.A + A(^ - 1) + huN ^''''^ 

2 4 ^2 ^12 

l-7 + 5/n27r C'(2) 1 1 1 ^, 1 ^ 

12 27r2 12 A 720A2 360A3 ^A^^' 

We add a few remarks. 
(i) The thermodynamic ansatz suggests that 

lnZN{p)--lnZNpil) = 0{N), (2.35) 

p 



while (trivially) Holder's inequality yields 

p lnZN{l) < InZNip). (2.36) 

(ii) One can also obtain a crude upper-bound as follows. From the symmetry of the 
integrand in Zn{p) ( p.l8|) we have 



Zn{p)=NI / l\i^^e-^^'^)\A{xo,...,XN.,)\'P. (2.37) 

J|xnl>lxil>--->lx«_il ,_n ^ 



'|xo|>|xi|>--->|xiv_i| ^^0 

In this domain 

N-l 



|A(xo,...,a:;v-i)|^^<2^^(^-i)n 



I |2p(Af-i-l) 



1=0 

The integral over the sector is smaller than the full integral, hence 

N-l 



Zn{p) < N\ 2P^(^-i) Yl {pk)\ , (2.38) 



k=i 



where the factor A^! 2^^^^ ^' is probably a gross overestimate. For further discussion 
of such bounds see the next section. 



3. Coulomb Gases. 

3.1. Quadrupolar mean field. 

In the interpretation of the norm (^]^) as a one component plasma canonical partition 
function, the gaussian measure in ( |2.19| ) is due to the mean potential created by the neu- 
tralizing background charge. This potential is not only determined by the locally constant 
density, but also by the shape of the bounding domain which, in the case studied so far, 
is a disk and therefore induces a rotational symmetry around the center, the minimum of 
the harmonic potential. Far from this minimum the global form of the potential has no 
consequence on the local correlations at a given density, neither on the residual thermody- 
namics. This suggests to add to the quadratic zz term an harmonic function (equivalently 
the real part of an entire function) which does not modify the constant charge density. To 
ensure convergence of the partition function we are limited to the real part of a quadratic 

9 



polynomial, hence up to a translation, to a multiple of a quadrupolar mean field V = Rez^, 
created by far away charges. 

It is rather remarkable that the plasma model with a quadrupolar component of the 
mean field can also be completely solved for P = 2 (s = 0). We give here a short description 
of the method and of the corresponding results. We want to evaluate the partition function 
depending on the intensity of the quadrupolar field that will be called th^ {n real) so that 
the total potential 

V{z,z) =zz- -th^{z^ + z'^) (3.1) 

confines the charges. To compute 

Z^{\-ii) = /" J]^e-^^^^"^ + ^*'^^(^'+^1)|A|2 (3.2) 

i 

we just need to be able to construct the polynomials Pn{z) {z G C) with real coefficients 
verifying the unitary orthogonal relations 

(fz 

■W{Z,Z) Pn{z)Pm{z) = Snm 

(3.3) 




TV 

w{z,z) = e-''+^"'^'^''+'"^ 
degiPn) = n = (0,1,2,...) 
lim Pr,{z) = (n!)-^/V 



We start from the identity 

" ^ -Z''AZ/2 



e-^ ^^/^ = = (3.4) 

TT ^y (a + d)'^ - 4cc 

where Z, Z' are the column and row vectors (^) , {z, z) and A = [ _ , J is a hermitian 
matrix such that {a + d)"^ > 4cc. Set 

zz--thix{z^ + z^) = -Z^AZ, (3.5) 

where 

10 



Now introduce the vector X = (^) and write the foUowing identity 

_ / ^ exp -i(Z - A-^xyA{Z - A-^X) 



=< exp (-{X^Z + Z^X) - -X^A-^X^ > 



(3.7) 



=< exp(x2; + xz 



1 - (tV)' 
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i-(thfiy 



Rescahng x as x — y — ^^^^ — 
on the left hand side, we get 
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Tj^u and moving the term exp( 
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TJ m 



(3.8) 



The r.h.s. of this expression involves generating functions for monic Hermite polyno- 



mials 



n 



exp(x'U- y) = ^Hn{x) — . 







Expanding both sides we have 



E' 
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thu' n\ ^^ l\k\ 

n=0 '^ k,l=0 



"'"'' <Ht{zJ^)H,(z 
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sh2n' 
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and since in the average < ... > the integral of the denominator is chfi, we get 

V2, 



file-''+-^"'''^'"+'"^Hk{zJ-^)Hi{z ' ^ 



sh2fx' 



The polynomials Pn are therefore given by 



sh2fx' 



chfi 



V. 



(thij) 



JkJ 



Pniz) 

The first few P's read 



(tM^/'^n(^) 



n!(cMi/2 

Po{z) = (cM"'/' 
Pi{z) = {chij)-^/^z 

F2(-2) = 1= — [z — shfx ch^j 



Psiz) 



+ 



-= — {z^ — 3shfj, chu z) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



11 



Substituting in the Vandermonde determinant \/ni{chn)^'^~^^''^' Pn{z) for z^ we readily 
evaluate the partition function as 

y /I n^ = (ch^) (3-14) 

the superimposed field is contributing only to the dominant term of the free energy by 
A^^ logchfj., this term being non extensive. After subtraction, the residual energy is there- 
fore independent of //, that is of the mean field at large distance as was to be expected. 

It is also interesting to compute the local density of the Coulomb system of N charges 
in that neutralizing background. We have 

N-l 



Pn{z, z) = W{Z, Z) ^ Pn{z)Pn{'i 



n=0 

N-l 



which for // = reduces to 



(3.15) 



z 



-^ cha n\ " Jsh2n' 

n=U ' 



|2 



^-1 / -\n 

PN{z,z)\^=o^e-^'J2-^- (3.16) 

^ — ^ n! 

n=0 

The asymptotic behaviour in A^ for \z'^\ of order A^ can easily be found outside a transition 
region, for which a finer analysis is required. Since for fixed z 

lim pis[{z,z) = l Wz, (3-17) 

N^oo 

there exists a curve Cn{p), being a continuous deformation of the circle \z\'^ = N, such 
that the asymptotic density pN is equal to 1 in the interior of the domain bounded by Cn 
and zero outside. To evaluate Cn asymptotically, we write that the variation of density 

Pn+i{z,z)- pn{z,z) (3.18) 

is maximal for z G Cn- 

Using the large A^ asymptotic behaviour of the Hermite polynomial Hn when ^ is 
finite 

^,„^M£l±^VE£.«,(2, + e- + l), (3.19) 

12 



with the definition of 



'"^'Tros^' ''''*>°'' ''■'°' 



the variation of density ( 3.18 ) is written, in the "exponential" approximation: 

(3.21) 



'zz + -thfi{z^ + z^) + N{lnthfi + 1 + Re{2(j) + e~^^)) 



exp 
The maximum over z is given by the vanishing of the derivatives with respect to z, z, i.e. 



-z + th^z + N^(l-e-^^)) = (3.22) 

oz 



and its complex conjugate. According to ( p.20| ), 



sh(f>^ = , -^ (3.23) 



so that (|3.22|) amounts to 

thix = \e-^^\. (3.24) 

Putting (f) = ^ + irj , the equation of the curve is 

z ^y/N{e^ cos r] + ie~'' sin r]) (3.25) 



which is an ellipse with semi-axis a = yNe^ , h = yNe~^. The area of this ellipse is 
TxN independently of /i the intensity of the quadrupolar field. The effect of an added 
quadrupolar field is therefore to deform the Fermi liquid disc into an ellipse with the same 
area. 

3.2. N charges on the sphere. 

The thermodynamic equivalence between the harmonic plasma and the system of 
charges on a sphere has been proven for (3 = 2 and is presumably valid for all temperatures. 
The spherical geometry avoids giving a special status to the center of the disk thereby 
increasing the symmetry and allows a clear definition of the pressure since the potential 
function does not depend explicitly on the density. 

The Coulomb potential created at point fl{9, (p) by a unit charge located at the north 
pole {9 = 0) of the unit sphere, is defined, up to a constant, by the Green function of the 
spherical Laplacian. The source term must be supplemented by a neutralizing background 
which can be taken to be uniformly distributed (to respect the spherical symmetry) and 
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the total charge of the source must vanish on the sphere, a surface without boundary. The 
potential satisfies therefore the inhomogeneous equation 

VuV = -2nd{u) + J (3.26) 



with 



f5{u)da = l (3.27) 



da = d{cose) A # , da = An (3.28) 

and the constant | represents the neutralizing density. Putting rj = cos9, we have 

B BV 1 

Up to a constant the solution is, 

V ^--ln^^ = -ln(sin-] =ln{-) (3.30) 

where r is the euclidian distance between the charge and the point on the sphere. We fix 
an additive constant requiring a zero average value of V{u), giving the potential 

2 1 f) 

Vo = ^^(-) - 7^ = -ln{\/esin-) (3.31) 

I Zi z 

which satisfies J V^da = 0. The correspondence between the spherical system and the 
original planar one is obtained by a stereographic projection from the south pole on the 
tangent plane to the north pole. The complex coordinate in this plane is 2z, with 

z = tg^-e'^ (3.32) 



da dz [\dz 1 (tz 1 



(3.33) 



47r 27rz {\^\z\^Y tt {\^\z\^Y 

while the angular distance 712 between two points z\ and z^ is given by 

2 -(l + |ziP)(l + |^2P) ^^-^^^ 

14 



On a sphere with unit area, the partition function at inverse temperature /J = 2p and with 
two-body potential ^(712) is defined by the average 






<^>- U'^^ 



^' ' (3.35) 

dai 



471 

For a sphere of area S, the partition function /^v is to be multiphed by 5""^. We can also 
write 

^-M-/n^(rTI-jIwn5Tini^.-%P''- (3.36) 



With the choice of potential Vo(7i2) defined in ( p.31|) , the canonical partition function 
is 

Q^(p)=et^(^-i)/^(p). (3.37) 

We note the exact results 

InW = TTn^ I n A (3-38) 



^ ^ \p=0 



and 



The partition functions on the plane Zn {p) eqn. ( ^^) , and on the sphere In (p) eqn. ( |3.36|) , 
differ only by the weight functions. They become equivalent in the central region \z\ = 

In appendix A, we find from (|3.37|) and ( |3.38|) the asymptotic free energy per unit 
charge F7v(l) (/3 = 2) 

_2F„(1, = 'JI^ = hn2.N - 2 + Oi'-f). (3.40) 

In contrast with InZ^ for the neutral plasma, we note the absence of terms in N'^lnN 
and A^^ in IhQn due to the choice of potential Vq having a spherical average zero. If we 
want to compare the two systems, planar and spherical, we have to renormalize the scale 
of the potential, which is equivalent to introduce an appropriate continuous background 
following the convention of Caillol and al. (formula 3.4 in their paper). 
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The residual free energy is now given by the following form 
-2/(1) = lim hn ^^^^^^ 



= lim ( -ln2TTN - -\ - - [InN - 1) (3.41) 

-log27r-lJ =-0.0810614... 

which agrees with the residual free energy of the neutral plasma i.e. the term proportional 

to A^ in (2.34). 



Remark on the canonical pressure 

The following section is devoted to obtaining bounds for j{p) when p > 1, to confirm 
the dependence in N observed for p = 1. On the sphere, it is the quantity 

g^(iv-i)/2 1 1^ |2y |A^,|2 = rrsin2 2^ (3.42) 

that produces the correct extensive property. 

Assuming for the time being that for any p the same InN term as in formula ( |3.40| ) , 
this gives the free energy per particle for A^ > 1 on a sphere of unit area, S" = 1. For S 
arbitrary we must add NlnS to InQ^. Although the energy of the system of A^ charges 
on the sphere is not extensive, a global pressure denoted by p not to be confused with the 
notation p, is defined by (3p = qs — ^^^"^ ~ ^ ~ density = p. The non-trivial term in 

( p.40|) is not contributing to the global pressure. On the other hand if we interpret the free 
energy Fis[{p) as that of an extensive system of density p, eq. ( |3.40|) gives a supplementary 
contribution to the pressure 

The resulting equation of state /3p = p{l — j) is that of a neutral plasma, a system that 
probably cannot remain homogeneous for /3 > 4. 
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3.3. Holder's and Hadamard's inequalities. 



As we noted in (|2.36|) , for the plasma on a plane Holder's inequality requires Fn{(3) to 
be a non-increasing function of /3, that is the "potential entropy" S{/3) = /?^^ is negative. 
This inequality reads 

< a" > > (< a >)^ Vs > 1 (3.44) 



where 



1 ^ 



(3.45) 



is an average value, (possibly a weighted average), of non- negative quantities. It follows 
that < a^ >i/^ is a non-decreasing function of /9 for (3 >0. We simply need to apply these 
relations to a = e~^°^, s = -^ > 1 



From this we derive the inequalities 



(3N 



lnZN{(3) non — increasing 



(3.46) 



F^m < FNi2), (3>2 
FNi(3)>FN{2), 0<(3<2. 

Similarly on the sphere, if the residual free energy is defined for all p by the limit 

-2f(v)- lim -In-^^^ 
we obtain 



(3.47) 



(3.48) 



(3.49) 



/(p) </(!), P>1 

/(p) >/(!), 0<p<l. 

Now let us look at Hadamard's inequality to obtain a lower limit to the free energy. 
We start from the identity 





N-l 


det 

i i 






n=0 



deti 



N-l 



EJV —J 
n=0 



CnZj Z^ 



n 



N-l 
n=0 '^n 



(3.50) 



where cq, ci, ..., cn-i are A^ positive indeterminates. The required inequality reads 



IA1 < 



2| CLi ■ ■ -aN-iGN 

CqCi ■ ■ -CN-l 
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(3.51) 



where the a^'s are the diagonal elements of the positive matrix, the determinant of which 
appears in the numerator of the r.h.s. of ( |3.5C1| ) 



N-l 
— \^ 

n=0 



2n 



z = l,2,...,iV 



(3.52) 



Define the polynomial of degree A^ — 1 



N-l 



(3.53) 



n=0 



Substituting in the expression ( p.36| ) for the integral In^p) the inequality ( |3.51| ) and writing 



Pj, we get 



In{p) < inf 



POO T-riV 
Jo lll=l 71 



«(Pj) 



dpj 



3 = 1 (l+p,)P(~-l) [l+PjY 



n 



N-l 
n=0 '^n 



that is 



In{p) < inf 



W (Rn CnY 



N 



(3.54) 



(3.55) 



with the convergent integral Ap{c) given by 



AJc) 



aP{p)dp 



(l + p)(^-i)?'+2 



(3.56) 



This integral is a homogeneous polynomial of total degree p in the N parameters [c]. 
Therefore , for every p > 0, we have 



InlNJp) 

N 



<inf 



inA,{c)-^J2^ 



incr. 



(3.57) 



The r.h.s. of (|3.55|) being homogeneous in [c] of degree zero, we have to seek the minimum 
according to the direction of the vector cq, ..., cn-i- In the first "quadrant" oi Rn {cn > 0), 
this continuous positive function tends to infinity on all coordinates planes. It therefore 
has an absolute minimum in an interior direction, and at that point it verifies 

1 dA pi 



; n = 0,l,...,A^-l 



That is 



A dcn N Cn 

"°° aP-^ipynP'^dp _ A 

(l + p)(^-i)P+2 ~N 
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, J., . . . , 



Vn. 



(3.58) 



(3.59) 



Prom homogeneity, only A^ — 1 equations are independent since the sum is given by Euler's 
identity. When p = 1, the solution is obvious since ( |3 .591 ) gives 

•" ""* --* (3.60) 



/„ (l+p)"+i N 
i.e. 

^n=A[^~^) (3-61) 

where the constant A is left undetermined. This yields the polynomial 

aip)=A{l + pf-^ (3.62) 

and we get the following upper bound 



- y In 



N - N ^^ ' V n 



n 

N-1 
i 
< - 



N "^ \ n 

n=0 



(3.63) 



or 



(nSn!)^ 



Let us denote by B]^{j>) the Hadamard upper bound (3.49) (given explicitly by the r.h.s. 
of (|3.64|) for p = 1). Comparing with the exact result ( |3.38| ) 



1—-^ 



from which we get 



(3.65) 



'^^'^^r^Oi'^) (3.66) 

This shows that the upper bound of —2Fn{1) given by Hadamard's inequality and the 
exact value are only differing by 1, meaning from (|3.41| ) 

/(I) = --ln2n +-> -hn27i. (3.67) 

Let us turn to the general case. The polynomial a(p) given in ( |3.62| ) for p = 1 is still a 
solution for arbitrary p up to a multiplicative constant. More precisely the ansatz 
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verifies the stationarity equations (|3.58|) . Indeed, after dividing by A 



dp= — (3.69) 



n 



[l + p)iN-i)p+2 '^ N 



'0 

Since tfie dependence on p drops out, tfiis is an identity in view of our previous discussion. 
The only doubt that remains is whether this stationary point is the absolute minimum. 
This is certainly the case in the neighbourhood of p = 1 by continuity if we note that 
the matrix Hn^m = CnCm-^ -^— Ap{c) is positive definite and that the function A^ [c) is 
convex in R^ hence admits only one minimum on the plane 'Ylm ^'^(^n = 0. 
The inequality ( p.57 ) can now be written 



^^^^^^<lnA-^yin{Ay^(^-^ 



N-1 ,^ ^. (3-70) 



~ N '^ \ n 

The upper bound for '^ ^^p' is according to ( |3.63| ) independent of p and the same property 
holds for 

Np Np 2 2 ^ ^ 

Consequently the residual free energy is bounded by 

/(I)- ^</(p) </(!), P>1. (3.72) 

where /(I) = —^ln27r + ^. Numerically we get 

-0.4594692... < f{p) < 0.0405307 (3.73) 

The lower bound is of no interest for < p < 1 since /(p) > /(I) in the high temperature 
region (i.e. p < 1). 

3.4- The zero temperature limit for finite N. 

Before going to the thermodynamic limit, i.e. keeping N finite, the free energy 
— ■klnlN(f3) of the system with A^ unit charges on the sphere has as a limiting value, 
when /? ^ oo (or equivalently p -^ oo), the absolute minimum of the potential function 

V^ = -^/n(sin^) (3.74) 
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or 

N 

2Vn = -/n|Ajvp + (A^ - 1) ^ ln{l + z^Zj) (3.75) 

With the definition (|3.48| ) of the residual free energy, we have for the finite system 

hm 2/(p) = 1 inf (2Va.) - ^^ + ^MN\. (3.76) 

This result remains valid in the thermodynamic limit provided the residual entropy van- 
ishes. Up to an overall rotation a stable equilibrium configuration does exist on the sphere 
since Vat reaches its minimum. The equilibrium configurations are solutions of the A^ 
equations ^ = (and c.c.) 

having in the planar limit the simplified form 

^ 1 
Y, —^, -{N- m = (3.78) 



■ , Zi- z 



From ( |3.77| ) we get the sum rule 



^ - AT 

ZnZn ]\ 



^TTTj-ri <^-™' 



j=i ■' ■' 

which reads in polar coordinates 

^ cos Oj = (3.80) 

j 

This result is true irrespective of the location of the pole, therefore we have ^ ■ Uj = {uj 

being the unit vector with polar angle 9j ) , which means that the center of gravity of charge 

is the center of the sphere as one would naturally expect. In the thermodynamic limit the 

relation (|3.79| ) allows to show that the radius of the confinement disk of the plasma is 1 

(with gaussian weight e"^^^^). 

The equations ( p.77|) are formally invariant under rotation (zj -^ ^j = _"f ^ ^- , aa + 

77 = 1). It is easy to find solutions for configurations invariant under a finite rotation 

group. As an example let us consider the case A^ = 12, where 2 configurations come to 

mind: 
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-12 (go) : The 12 charges are at the vertices of an icosahedron. The order of the 
symmetry group is 60. 

-12(24) : The 12 charges are at the vertices of a cuboctahedron (center of the edges of 
a cube). The order of the symmetry group is 24. 

From symmetry we see that the electric field created by the 11 charges on the twelfth 
located at the pole is zero: '^^ j- = and the equations ( |3.77|) are all satisfied. Turning 
to the question of stability one finds that the configuration with highest symmetry group 
has lowest energy. In appendix B we show that the configuration of the cuboctahedron 
is unstable. By a similar method we can show that the icosahedron is stable, that is the 
Hessian is semi-positive. The Hessian cannot be strictly positive since it admits three zero 
eigenvalues corresponding to the three generators of global rotations. The equation for the 
vibration modes A, (A = u"^) is written as 

E /'!% +(N- l)(^"^(l -^)- ^'^0 - (3.81) 

J \Zi Zjj 

for the 2N eigenvalues {5i,Si). We exhibit in appendix B one ternary unstable mode of 
the cuboctahedron. 

Let us briefiy present the method of calculation of the potential function Vat for some 
symmetric polyhedral configurations belonging to the octahedral group or the icosahedral 
group with N = 8, 12, 14, 20, 30, 42. We will also give a table of residual energies to compare 
with the previous bounds. 

Given formula ( p.74|) we only need to know distances between the vertices of the 



polyhedron, only 3 distinct ones in the case of the icosahedron. For higher configurations 
it is more appropriate to compute systematically the discriminant A^ (formula (|3.75|) ) 
knowing the invariant polynomials (having roots Zj in the complex plane corresponding to 
the vertices on the sphere). Let us take for instance the configuration 12(24) of mid-edges 
of a cube (or the cuboctahedron) . Klein |^ gives the following polynomial 

f{z) = z^^ + l- 33(zS + z^) 

(3.82) 

with 

a^ + a~^ = 6, a~^ — a = 2 
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We get a^ = 3 — v^ and a ^ = 3 + v^. We have A^ = Ylj^i f'{zj)- The roots appear in 
four quartets. We give the value /' for a typical root, then the product over the quartet 

n /'(e^^/^) = (4 X 6^)4 

4terms 

/'(a) = ^a^{a^ + l)(a^ - oT^) = 4a^ x Ga^ x 6 x -2v^ 
J] /'(a) = (a^ 23v^ X 6^)4 

4ierTns 

52\4 



JJ /'(a-1) = (a-^ 23v^x62 



also 



Aterms 

A2 = (4 X 6^ X 2^ X 2^ X 6^)^ = (2^^ 3^)^ 

— /n|A|2 = 6.125058601 
12 ' ' 



(A^ - 1) ^ ln{l + \zj\^) = 11 X 4(/n2 + /n(4 + V8)(4 - ^8)) 
i 

= 11 X 2^/n2 

— {N-1)Y^ ln{l + l^jf ) = 10.16615865... 



Ap (217 36)4 



„-2Vi2 _ 1^1 _ V^ "* ^ _ /9-9q2Nl2 

^ ~ TT /I , I., I5IM1 ~ olfivll ~ '^^ "^ ^ 



^ _i_ ^ |2U1 916x11 

2Vi2 



-2/n3 + 9/n2 = 4.041100047 



12 

The same method applies in the following cases: 
-12(60)- Vertices of an icosahedron 



with 



f{z) = ziz'"" + 11^^ - 1) = ziz'' - a^)(^^ + a"^) (3.83) 



a-^ -a'' = 11 

27r 

a = 2 cos — = 2^ — 1 

5 

a-i = 2 cos- =2^ 
5 

I A|2 _ 1-25 



2Vi2 = 55/?i5 - 25/n5 = 30/n5. 
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-20(60)- Vertices of the dodecahedron 



f{z) =z^^ + l- 228{z^^ - z^) + 494z^^ = {z^ + a^){z^ - a-^)(^^ + 6^)(z^ - 6"^). (3.84) 

with 

a~^ -a^ = 114 + 50^5 

6-5 -b^ = 114 - 50^5 

< a < 6 < 1, 

1^12^2603^0595 

2Vi2 = 190/n(a + a~^){b + b~^) - (60/?i2 + 10/?i3 + 95/n5). 
-30(60)- Mid-edges of the dodecahedron (or the icosahedron) . 



f{z) = ^30^;^^ 522(^25 _ z^) - 10005(^20 + ^10^ 

= (^1° + l){z^ - a^){z^ + a-^)(^^ + b^){z^ - b-^) 



(3.85) 



with 



a~^ -a^ = 261 + 5^^5 
6-5 _ 6^ = -261 + 5^v^ 

I A |2 _ 290o60r205 

2Vi2 = 290/n2(a + a"^)(6 + 6"^) - (90/?i2 + 60/?i3 + 205/?i5) 
-42(60): Configuration with the 12 vertices and the 30 mid-edges of the icosahedron. 

|A|2 = 29036°5305 

For this mixed configuration the discriminant is obtained by multiplying |Ai2p x IA30P 
by the factor Yliif fsoi^i), the product being on the 11 roots of /12 which gives us the 
extra factor 

(11^ + 4)^(261 + b^Vb - 11)^°(-261 + b^Vb + ll)^^ = 5^^ 
2^12 = 410/n2(a + a-^)(6 + 6"^) - (lOO/nS + 60/n3 + 90/n2) 

The results for the free energy per particle ^ and the residual free energy 2/ (00) = 
^ — -^^^ + -^InNl are shown in the following table. We note that the approximate value 
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/(oo) ~ -0.32 is in the limits defined in ( p73|) namely : -0.46 < -0.32 < 0.04. Finally 
Holder's inequality restricts the width of the interval from 0.5 to approximately 0.37. 



A^ 


A^-i/n A|2 


X 


F = N-^V 


Y 


-/ 


8(24) 






1.1280580 


1.4186061 


0.290548 


12(24) 


6.1250586 


10.166158 


2.0205500 


2.3335996 


0.313049 


12(60) 


3.3529956 


7.376590 


2.0117973 


2.3335996 


0.3218023 


20(60) 


10.273577 


18.081646 


3.904034 


4.2208047 


0.3167702 


30(60) 


15.274491 


27.880218 


6.3028631 


6.627848 


0.3249848 


42(60) 


14.742346 


33.202164 


9.2299092 


9.548976 


0.3190676 



Table I: Residual free energy at zero temperature on the 
sphere for some symmetrical configurations. 

( with X = {1- A^-i)E/n(l + \z\^) and Y = \{N - 1 - N-HnN\)) 



As already mentionned, we show in appendix B that for A^ = 12 the configuration of 
the cuboctahedron is unstable. Its potential energy is larger by an amount 0.008752... than 
the one of the stable icosahedron configuration. This energy difference can be interpreted 
as an "activation energy" associated to the cuboctahedral configuration as compared to 
the icosahedron. It is interesting to note that the associated inverse temperature Pc = 
2pc = 1/0.0087... ^ 114 is of the order of magnitude of the estimated fusion temperature 
corresponding to the transition from a quasi-solid phase to a disordered phase. 



4. Expansion in characters 

We return to the expansion of Laughlin's wave functions in terms of Slater determi- 
nants. Using the notations of section 2 we would like to get a better understanding of the 

(s) 1 

coefficients gy ■ For short, when s = 1 or p = 3 corresponding to filling fraction ^ we set 

(1) - 
9y = 9y- 

Consider the monic polynomial of degree N, P{x) with roots x^, < z < A^ — 1, 



N-l 



N 



p(x) = n ( 






Et-D' 



akx^ ^ ,ao = l. 



(4.1) 



i=0 



k=0 
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It has a double root if and only if its discriminant vanishes or equivalently if it has a 
common root with its derivative 



N-l 



P'{x)= Y,{-lfak{N~k)x 



N-k-l 



(4.2) 



fc=0 



The discriminant, as a symmetric function of the roots is obtained by eliminating x between 
( [4.1| ) and ( |4.2D . A systematic way to perform this elimination is to think of the powers of 
X as independent variables and to consider the polynomials, 



x'^Pix), 0<m<N -2, 
x^'P'ix), 0<n<N-l, 



(4.3) 



as 2N — 1 linear forms in the 2N — 1 variables x^^ x^, ■ ■ ■ , a;^^~^. A root common to P{x) 
and P'{x) then requires the vanishing of the corresponding determinant which is therefore 
proportional to the discriminant, in fact with our conventions equal up to a sign. For 
instance for the classical cases, 



P{x) = X^ — 0"iX + 0"2 

D 



1 


-CTl 


0-2 


2 


-CTl 








2 


—a 



4<T2 - CTl 



P{x) = x^ — aix"^ -\- a2X — a^ 

1 -CTl 0-2 -(73 

1 —CTl a2 —(y^ 

D = 3 -2ai a2 

3 -2ai a2 

3 -2ai a2 

-afal - 18cricr2cr3 + Aafas + 4(j| + 270-3 



(4.4) 



(4.5) 



and in general. 



D = 



-0-1 
1 



0-2 
-ai 



N -{N-l)ai {N-2)a2 
A^ -(A-l)ai 



(-l)^aA. 

(-1)^-Va._i (-l)^a^ 



(-1) 





Af-l 



ajv-1 



N 



-(A-l)ai 



(-1)^-V^_i 
(4.6) 
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Alternatively knowing that 

N-l 
P'{x) = y^ {x — Xq) ■ ■ ■ {x — Xi) ■ ■ ■ {x — Xtv-i) 



i=0 



we have, P'{xi) = Hj.j^iixi - Xj), and, 



N-l 



J^i 



0<i< N -1, 



d{xo:---,xn-i) = Yi n "^^^ ~^^'^ ^ n ^'( 

0<i<N-l j;j^i i=0 

In particular for A^-th roots of unity Xi = exp{ ^ j^ — - 

P{x) =x^ -1 

N-l 

D{xo, ■ ■ ■ , X;v-l) = iV^ n (^^)'^"' = (-l)^-'A^^ 

i=0 

One verifies that this does agree up to a sign with the above determinant. 
From equation ( |2.6|) we have also 



iV(iV-l) 



1 1 

xq Xi 



^N-l ^N-1 

Xn Xi 



XN-i) = A^(a;o,- --^xn-i] 



1 

XN-1 



X 



N-l 
N-l 



X 



1 Xo 

1 Xi 

1 Xn- 



X 



,N-1 
%-l 



X 



N-l 
N-l 



(4.7) 



(4.^ 



(4.9) 



(4.10) 



where, 



to ti 

tN-1 ■■■ 



tN-1 

In 

t2N-2 



tk — Xq ~r X^ + ■ ■ • + Xj^_-^ 



P'jx) 
P(x) 



E 



tk 



fc=0,|x|>|xi| 



X 



k+1 



(4.11) 



In particular to = A^. For a generalization (due to Shiota) of formula (|4.10| ) for all s see 

appendix C. In the example of A^-th roots of unity tk = except when /c is a multiple of 

A^, in which case it is equal to A^, the above reduces to 

10 ■•• 

1 

N 1 



Af(Af-l) 



(-1) — 5 — D^N 



1 



AT (N-l)(N-2) 

N^i-iy — ^ — - 



(4.12) 
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in agreement with (O). 

Let us apply ( [4.1UD to give a concise - but not very effective - formula for the coef- 
ficients qy in the expansion over Schur functions. For this purpose we recall Frobenius 
reciprocity formula expressing a product of t's as a combination of Schur functions chy 
having coefficients xy, the characters of the symmetric group of order |y| pertaining to 
the same set of Young tableaux. Recall that xy is only a class function on S\y\- We denote 
the classes (1)"S (2)"^, ■ ■ ■ , i.e. cti cycles of length 1, 012 cycles of length 2 , etc... then 



j.ai .02 
'■I ''2 



J2 XY{{in'2r---)chY{xo,---,XN-i) 

Y:\Y\=Y,kak 



(4.13) 



For each class (l)"i(2)"2 ... we have a number - in fact a relative integer - that we write 
Xy((l)"^(2)"^ ■ ■ ■). Define symbols (1)(2) ■ ■ ■ that can be multiplied and linearly combined 
(with coefficients in 2 , Q ,1R or (D). Attribute to (k) the degree k and consider any linear 
combination of fixed degree |y|. Then we can extend the definition of xy by linearity 



ai+2a2 + -- 



"01,02, 



.(1 



i«l/ 



\Ct2 



J]a«,,„,,..Xy((l)"n2; 



02 



(4.14) 



\Y\ 



If the coefficients are integers so will be the value of xy- With this notation and for s = 1, 
Qy, we find by combining ( |4.1UD and ( |4.13[ ), 



9? 



A2(xo, 



XAr_i) 



\Y\=N{N-1) 



N 
(1) 



(1) (2) 

(2) (3) 



(iV-1) 



(N-l) 
(AT) 

(2A - 2) 



chY{xo,- ■■,xn-i), 



(4.15) 



where the sum on the r.h.s. extends over Young tableaux with N{N — 1) boxes, at most 
N lines (otherwise c/iy(xo, ..., X7v_i) vanishes) and at most 2{N — 1) columns (since the 



leading term in A'^(xo, 



xn-i) is A(x|], 



' ^N-l, 



We call these Young tableaux 



weakly admissible and generalize this notion for arbitrary s > 1 to those tableaux of area 
sN{N — 1) contained in a rectangle of size 2s{N — 1) x N. The reason for the qualifier 
weakly is that some of the tableaux have necessarily vanishing coefficients as will be shown 
below. In any case for (weakly) admissible tableaux, 



9y = Xy 



N 
(1) 



(1) 
(2) 



(N-l) (N) 



(N-l) 
(A) 

(2A-2) 



(4.16) 
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Even though this formula solves in principle the problem of finding an explicit expression 
in the case s = 1 it is not very manageable. Indeed the r.h.s. is a linear combination of 
A^! values of characters of >S'7v(Ar-i), a formidable task to evaluate. It is however possible 
to dispense altogether with the characters xy of the symmetric group as follows. The 
characters chy of the linear group are expressed in terms of the normalized traces 9k = ^ 
(gJ|) as 



ch 



Y{fo,...,fN) 



Efc>iA;afe = 



PfoiO. 

xviiirH'^r' 



(4.17) 



\Y\ 






Here Pn(^i) stands for the elementary Schur functions, the characters of the n-th symmetric 
power given by 

nai f\a2 

PniO.)= E -h-h--- (4-18) 



The * in the determinant indicate that moving a unit step to the right (left) the index of p 
increases (decreases) by unity, and one agrees that po = 1 as well as p„ = for n negative. 
The second part of ( |4.17| ) (from which ( [4.18| ) follows) is obtained by inverting eq. ( |4.13| ) 
using the orthogonality and reality of the Xy 's. For short set dk = -^ in such a way that 



dkPnid.) =Pn-k{0.) 



(4.19) 



Using these notations we rewrite ( [4.16| ) as 



9Y 



d 



N di 
di 82 

'n-1 Qn 



On-i 




On 






X 


%N-2 





Pf.^Ad-) 



Pfoid- 



(4.20) 



This gives a more manageable way of computing these coefficients. Notice that here we 
treat the ^'s as independent variables and that the amount of derivatives is such as the 
right hand side is a pure number (indeed an integer). 
Finally if we observe that 



exp^z^OkPniO.) 



1 f d 



n\ V dz 
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exp 2^ z Ok 
1 



(4.21) 



we can even dispense with the p^s and reach a simple expression 



gr = 



d. 



fjv_i-JV+i 



d. 



lpf-l-N + 2 



-N+iy. (iN-i-N+2y. 



d. 



lO-N + l 



d. 



lO-N+2 



{lo-N+l)\ 



(lo-N+2y. 






X 



to 



In-i 



N-l 



t 



2N-2 



(4.22) 



As before to = N and on the right hand side the terms -^d^ where k < or k > 2N — 2 
are set equal to zero. 

From these expressions we derive a number of properties, starting from the most 
obvious ones which follow immediately from the definition, with p = 2s + 1 
Property 
(i) All coefficients gy are integers. 

(») sit,.-. - "'" 



9p{N-l)-lN--l,p{N-l)-lN-2,---p(N-l)-lo 

^ (s),Af 
lN=pN ~ 9ig^...iN_:^ ' 



.■■■) is),N+l {s),N 

(ill) gl'' - "' " 



yo,io+p,---iN-i+p y^ov-'iv-i 



Assume that s = 1 and suppose that Y is contained in an iV x A^ square (Its area 
being N{N — 1) this is indeed possible). This requires fN-i < N or In-i < 2A — 1. If this 
is the case Y and its symmetric Y with rows and columns interchanged, are both weakly 
admissible. Since xy and Xy differ only by the sign of the permutation and since a cycle 
(k) has parity (—1)^, we have. 



9y = Xy 



N 
(1) 



;i) -(2) 

■(2) (3) 



(-1)^-2(A-1) 
If we set with A and B integers, 

then, 



(-l)^^-^(A-l) 



-(2A-2) 



gy = AN + B 



(-iy^g^ = -AN + B 
where the extra index, N or N+1 refers to the number of variables. 



(4.23) 



(4.24) 



(4.25) 
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consequently gy — (— 1)^(7^ = 2 AN and, 
Property 1 

For Y and Y weakly admissible gy = {—^)^9y mod 2N. 

Example: 

A (xo^xi) = chojixo^xi) — 3chr\{xo,xi) (4.26) 

and 1 = —3 mod 4. For more convincing cases one may look at the tables in appendix D. 
Now we have a slight refinement of Property 0. Suppose that p = 2s + 1 is an odd 
prime, then from Fermat's theorem (x^ — Xj)^ = x^ — x^ mod p. Thus, 

AP(xo,---,XAr_i) = A(x^,---,x^_J modp (4.27) 

and 
Property 2 

Apart from g^^ ^ (n-i) ~ ^' ^^^ coefficients gi ... ;^_^ vanish mod p if p = 2s + 1 
is prime. 

That g^ (n-i) ~ ^ follows from the definition (|2.10|) there being a unique term 

XoX^---x^_j^ with coefficient 1 in the expansion of Y[i>ji^i ~ ^j) ■ -'-^ particular 
Property 2 applies for s = 1 where apart from the leading one, all coefficients are multiples 
of 3 (see appendix D). One may generalize Property 2 when p is composite using that if 
P = P1P2 with pi prime then A(xo, ■ ■ ■ , xn-iY = A(xq\ ■ ■ ■ , x^_^)p^ mod pi. 

The expansion ( 2.1U ) over Slater determinants ranges from the most extended one 



with the sequence of /'s being (0,p, 2p, ■ ■ ■ , (A^ — l)p) to the most compact one leaving a 
huge hole around the origin (and physically very unlikely) with every level from s{N — 1) 
to (s + 1)(A^ — 1) occupied i.e. the sequence 

lo = s{N - 1), /i = s(iV _ 1) + 1, . . . , l^_, = [s + l){N- 1) (4.28) 

For general s let us describe the admissible sequences of /'s i.e. those for which the corre- 

(s) 

sponding coefficient gi ... i _ is not a priori zero. Subtracting successively 0, 1, ■ ■ ■ , A^ — 1 
we get the admissible sequence of /o < /i < ■ ■ • < fN-i of Young tableaux. 
Property 3 

For s > 1 the set of admissible Young tableaux is of the form: 
fk = 2sk + rik+i -Uk, < k < N -1 
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or equivalently, 

Ik = (2s + l)/c + Uk+i - Uk 

with non negative integers rik, such that 

no = UN = and < Uk < ^{uk+i + Uk-i) + s. 

The last inequality follows from the requirement that fk > fk-i- Clearly all admissible 
Y are weakly admissible since fo<---< /n-i < 2s{N — 1) — un-i and the last bound is 
smaller or equal to 2s {N ~ 1). 

To prove this property we consider for odd p 

A(a:o, ■ ■ ■ , XM-iT = U^x. - XjT = U(^~^r^ x^^-^^x^^^"') ■ ■ ■ x°. (4.29) 

For i > j write — = -^^^^^ x -^^^^ ■ ■ ■ — ^^ and expand the last product in ([4.291). One 



obtains a sum of terms of the form 



XN-2 .nM^i( ^N-3 .nM-2 ( ^0 .^ ^ p{N-l) piN-2) Q 

{ ) { ) '''\^ ) ^ -^N-l •^N-2 ■■■•^0 

Xn-i Xn-2 Xi 

— V{N -^)-'n-N -1 p{N -2)+nN -i-riN -2 pk+rik + i-rik m 

— ■^N-1 •^N-2 '''-^k ■■■•^0 



(4.30) 



with integer coefficients and non-negative powers. Both sides in ( [4.29| ) being antisymmet- 
ric we will find on the right hand side for each monomial the antisymmetric sum of its 
permutations. To identify these determinants it is sufficient to retain those monomials for 
which the exponents pk -\- rik+i — rik form a strictly increasing sequence. This is Property 
3. We denote by A^^ the number of admissible tableaux. 

For 2 < A^ < 29 we present in the following table the total number A)^ = An of 
admissible tableaux for p = 3, s = 1. 
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N 


Number A^ 


Ratio An/An- 


-1 


2 


2 


2.00000 




3 


5 


2.50000 




4 


16 


3.20000 




5 


59 


3.68750 




6 


247 


4.18644 




7 


1111 


4.49797 




8 


5302 


4.77227 




9 


26376 


4.97472 




10 


135670 


5.1436 




11 


716542 


5.28150 




12 


3868142 


5.39834 




13 


21265884 


5.49769 




14 


118741369 


5.58365 




15 


671906876 


5.65857 




16 


3846342253 


5.72451 




17 


22243294360 


5.78297 




18 


129793088770 


5.83515 




19 


763444949789 


5.88201 




20 


4522896682789 


5.92432 




21 


26968749517543 


5.96271 




22 


161750625450884 


5.99770 




23 


975311942386969 


6.02972 




24 


5909549998347426 


6.05913 




25 


35966989049703188 


6.08624 




26 


219805620101524672 


6.11131 




27 


1348411202206287872 


6.13456 




28 


8301060310442080256 


6.15617 




29 


51270095697195679744 


6.17633 





Table II: Number A^ of admissible tableaux for s=l. 



The last column gives the ratio of successive terms. The above reasoning does not 
however insure that this is exactly the total number of terms in the expansion of A^ in 
characters as some coefficients might still vanish. However experience up to N=5 (appendix 
D) seems to indicate that these accidents do not happen. It is curious to notice that the 
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condition Uk < ^{uk+i + Uk-i) + s is a discrete generalization of subharmonic functions. 
The number of admissible tableaux is smaller or equal to the number of weakly ad- 
missible ones 

^N,2siN-l){sN{N-l)) 

where n7v,M(/) = ^M,N{f) stands for the number of partitions of the integer / in at 
most A^ (non-empty) parts of size at most equal to M. Equivalently this is the number of 
tableaux in a M x A^ rectangle of area f. It is given in terms of the generating function, 

E^^-Mf) -^^ = (1 -,)(!_, ,)(!-, ,.)... (!-.,«) ("D 

Recall that in terms of an infinite sequence of variables 0i,92, ■ ■ ■ the (elementary) Schur 
polynomials (in the ^'s) are defined through, 

oo oo 

exp J2 ^"^n = Y, ^>- (^0 (4-32) 



Since 



oo 



we find, 



-^ = exp yx"" -^ (4.33) 

1 1 _ „A;(M+1) 

Y, ^nMIW = PNiOk = ^ ^ _ , ) (4.34) 

/>o ^ 

An elementary fermion-boson equivalence gives an alternative generating function for the 

restricted number of partitions, namely 

Y nA.,M(/)y V^'''^'^ = {l + y)il + yq)---{l + 2/?^+^"^) (4.35) 

N,f>0 

Hence 

Y^nMDq^^^-^ =PN{Ok = ^-^ /_ ^, ) (4.36) 

thus 
Property 4 

The number of admissible tableaux is smaller than the number of weakly admissible 
ones given by Il]\[^2s{N-i){sN{N — 1)). 

34 



Property 3 suggests to write the expansion over characters using the operator 

A(t,)= H (l-r,r,+i---r,) (4.37) 

l<i<j<N-l 

where the "box shifting" operatoiS n acts on characters as 

Ti chi^^,...,iN-i = chi^,...,h-i+i,h-i,...,iN-i- (4.38) 

The expansion of A^** over the characters is also obtained by retaining the admissible part 
of the action of A^*+^(ri) on the leading character c/io,(2s+i),2(2s+i),...,(2s+i)(Ar-i)5 namely 
that of the monomials r]^^...rjY^^^ which produce a strictly increasing sequence of indices 

Til < (2s + l)-ni + n2< ... < (2s + l)^-^ - nN-2 + riN-i < (2s + 1)^-^ - n^v-i, i.e. 



2s+l 



in) 



adm 



"!,■■■, "JV-l 
admissible 



cho^(2s+l),...,(2s+l}{N-l) 

E^(s) r m njv-i 7 1 

'-^{ni,...,njv-i} Fl ■■■^Af-1 C/^0,(2s+l),...(2s+l)(^r-l)J ^439^ 

Y 

where the n's in the summation are as in Property 3, and C\^\ = Qy if Y is obtained 
from the tableau {0, (2s + 1), ..., (2s + 1)(A^ — 1)} by action of a monomial t"^^ ■■■t^^Si ■ 

les us 
pansion of /S.{TiY^^^: 



is) 

This provides us with another compact formula for the ^r^ , by straightforward ex- 



l<i<j<iV-l 



{n} - / . \ -I I I \ .__ / I I \ „_ V- .__ V- „__ / ' 



iV-1 

(4.40) 



where the sum over rriij = rriji is restricted to the values for which the combinatorial 
factors make sense (( ^^ ),for < n < 2s + 1) we list here the first two cases 

N=2: 

'r^=Ci^^ = {-ir(^"^% 0<n<s (4.41) 



9y' .. . . , 

\ n 



The elementary action of Ti on a Young tableau is to move a box from line N — i to line 
N — i -\- 1. The actions of various r's commute with each other. 
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N=3: 



with respectively Y = {n, 2s + 1 — n} , < n < s for N=2 and Y = {ni, 2s + 1 — tii + 
77-2, 4s + 2 — 77,2}, 0<ni,n2; 2ni— n2<2s, 2n2 — ni<2s. More coefficients are presented 
in appendix D. The expression ( [4.4C1| ), although increasingly complicated with A^, provides 

(s) 

us with a factorization property of the gy 
Property 5 

For gy — C*!^ ^ -, as above, we have the factorization property: 

Cf ^ 1 = Cf ^ 1 X Cf ^ T (4-43) 

{no,.--,n-j-l,0,nj+i,...,njv_i} {no,...,nj_i} {nj + i ,...,njv-i} ^ ' 



This is a straightforward consequence of eqn ( |4.40| ), as the vanishing of Uj implies 



that of all niki, k < j < I, therefore the sum over rriij > breaks up into a product of two 
sums over niki with respectively k < I < j and j < k < I. 

Let us derive for p = 3 yet an alternative expression for the coefficients gy starting 
from equation ( [4 .81) which can be rewritten 



{-1)^^^^^A^{xq, ■■■, XN-i) = |P'(x)x°, P'{x)x\ ■■■, P'{x)x^-^\ (4.44) 

Recall that this is understood as a determinant with xq substituted in the first line, xi in 
the second, ... Introduce an other set of variables //q, ■ ■ ■ , fJ^N-i- One readily sees that 



iV(iV-l) 



(-1)^^ A^(xo, ■ ■ ■ , XN-l) = 



a^ ,„, . n „, . .._., (4-45) 



duodni ■ ■ ■ diiN- 



tj-k=xk 



-m, WN 



since xq, ■ • ■ , a^iv-i are the roots of P. Set for simplicity (prn = (~1) ^ctn 
o"o = 1) in such a way that 

P{^)= E (-l)'^fc^'^"' = E '^-^"^ (4.46) 

0<A;<Af 0<m<Af 
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iV(Af-l) 



{-l)^^A%xor--.XN-i] 



d 



N 



diiQ- --duN-i 



E 



bmo(l)m,-l ■ ■■4>m^_,-N+l I/"""" " " " /U"^^"^| 



mor'-yrriN-i 



t^k=Xk 



d 



N 



diiQ- --duN-i 



E 



m()<---<mjv-i 



'mi 



^mi — 1 



''Tnjv-i Y^Tnjv_i — 1 



-'mo-iV+l 
^mi-Af+l 



'miv-i-Af+1 



X |^""o...^"'~-i 



l^k=Xk 



^ mo...mAr_i 

m,o<---<"T.jv-l 



■^mo *r mo — 1 



''mjv-i V^mjv-i — 1 



We now change ruk — ^ rn/c + 1, substitute 
with the result that 



^mo-N+l 



^ruN-i-N+l 



V I^Tno-1 mjv-i-l| 



(4.47) 



(— 1) ^aN-rm and coUect the signs, 



A3(xo, ■ ■ ■ ,XN-i) = Yl i-^f'^Hmo + 1) ■ ■ ■ {mM-1 + I] 

mo<---<m,]v_i 
O'N-l-mo O-N-mo ■ (^2{N-l)-mQ 



X 



(^N-l- 



■mjv_i 



C"2(Af-l)-mjv_i 



I^TTio ^miv_i| 



(4.48) 



Set 






(4.49) 



Define the Young tableaux 



Y{f) = {0<fo<fi<---<fN<N-l} 
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(4.50) 



and Y{f") given by the sequence < /q' < ... < f']^_^ <N-1 with 

f^_^=N-l~mo=N-l-fo 



f'j(j.k = N-l + k-mk = N-l-fk 



(4.51) 



f(^ = 2{N-l)-mN-i=N-l-l 



N-l- 



Then 



U.'Tio . . . j."ijv-i 



A(xo,---,XAr_i) 

while according to a standard property of characters 



chY{f){xo,- ■ ■ ,xn-i) 



(4.52) 



O'N-l-mo 



(^2(N-1) 



-mo 



aN-1- 



■mjv_i 



'^2(Ar-l)- 



■mjv_i 



CTf" (Tfii 11 

JJV-I J]V-l^-^ 

(J r;/ _-|^ 



^/''-(7V-1) 






^/o' 



cV77^)(^o, -,a;7v-i^ 



(4.53) 



where Y —>■ Y corresponds to the transposition (i.e. the interchange of hnes and columns) 
of Young tableaux (a notation already used) in which elementary Schur polynomials {pk = 
trace of /c-symmetric power) get interchanged with elementary symmetric functions (a^ = 
trace of /c-antisymmetric power). This is recorded in the following picture where Y{f) and 
Y^{f') = Y{f") are two complementary parts of an A^ x (A^ — 1) rectangle. 

Fig. 1: Notations used for Young tableaux 



This gives the final result 

YeNxiN-1) 



JV(JV-l) 



+E/- 



(4.54) 



X (/o + l)---(/Ar-l + A^) chY{xo,---,XN-l) chy-ixo, ■ ■ ■ , Xn-i) 
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where the short hand notation Y ^ N x (N — 1) indicates that the tableau is inscribed in 
the rectangle A^x (A^ — 1) and if Y is given by the sequence of < /i < • • ■ < /n-i < A^ — 1, 
Y^ corresponds to the complementary sequence /q = </(<•■■ < /jv-i ^ -^• 

As a (trivial) example using the notation < /q, /i ■ ■ ■ , In-i > for c/i/g=/Q^... j^_j=z^_^_7V+i 
we have, 

A^(xo, xi) = —2chjcD + 3c/iD c/iq — 6c/in 

= —2ch\jj + 3{ch\jj + chn) — Gchn 

= chm - 3c/in =< 0, 3 > -3 < 1, 2 > 
and the slightly less obvious one, 

A^(xo, xi, X2) = —Qchrrn + 8cho c/j- rrn — 10ch\jj chrn 
—12chr\ ch rrn + Ibchrn chrn — 20chrn chn + 24chn ch\ 
— 30c/i rn ch\jj + 40c/im ch\j — GOchrn 



rm 



-6c/irm + 8(ch rrrn + ch rrn + chrrn) - 10(c/irrm + chrrn + ch 



- 12 (c/l rrm + ch rrm + c/j- rrn ) + ISCc/l rrm + c/j- rrm + c/j.|-pp| + 2c/lmi| 



-20(c/irm + chrrn) + 24c/j. r-n-n - 30(c/i n-rn + c/irm) + 40(c/irrp + c/im) 



-60c/ir 



=< 0, 3, 6 > -3 < 1, 2, 6 > -3 < 0, 4, 5 > +6 < 1, 3, 5 > -15 < 2, 3, 4 > 
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Since for given A^ the Clebsch Gordan rule yields, 



Chy, Chy, = J2 ^vlv^chy, (4.55) 

V3 



til J V 

as 

2 



with Ny^ y non negative multiplicities, and II3I = |Yi| + 1^21, we rewrite equation ( [4.54| ) 
A^(xo,- ■■ ,xn-i) = 



N(N^l) I ^JV-1 



E E Nl,Y,i-^)'''^^''^^''''^'\foiY^) + l)■■■ifN-^iY^)+N) 

Y \YieNx{N-l) 
X chy{xo, ■■ ■,XN-l) 

(4.56) 
Equivalently, 

gy= Yl Nl^y.{-l)''^'^+^o''f^(''^\f,(Y,) + l)...{f^_^iY^)+N) (4.57) 

YiENxiN-l) 

The computation of gy is here reduced to the application of the Hall-Richardson rule 
to yield the Clebsch Gordan series over the Af = J2f^N,N-i{f) products of characters 
pertaining respectively to a Young tableau inscribed in an A^ x (A^ — 1) rectangle and its 
complementary. 



5. Sum rules 

A famous identity due to Dixon (1891) states that P 

A considerably generalized form first conjectured by Dyson was proved by Gunson and 
Wilson independently. Let as before, 

P{x)= J] (x-x,) (5.2) 







0<i<N-l 






Set 












V{xQ...XN-i;ao...aN-i) = 


y. (P'{X,)Y 


= 11 (1- 

0<i/j<Af-l 


^ J \ a. 




11 I ^N-1 ) 
0<i<N-l ^ -^i ^ 



(5.3) 
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Consider the constant term in V understood as a Laurent series in the variables Xi 

G'A.(ao,---,a^_i)= / W ^P(e^^°,...,e^^--^ao,...,aAr_i) (5.4) 

•^ 0<i<N-l 

then 

Theoreni(Dyson, Gunson, Wilson ): 

G'Ar(ao, ■■ ■,ajv-i) = „ V (5.5) 

Dixon's identity is nothing else than the particular case where A^ = 3 and ao = ai = a2 = p- 
The remarkable proof due to Good Q goes as follows. Remarking from Cauchy's formula 
that 

we multiply V{xo, ...xn-i', ciq, ■■■(in-i) by 1 written as above to get for the constant term, 

G'Ar(ao,- ■ ■,aAr_i) = ^ G'Ar(ao, ■ ■■,ai - 1,- ■ ■,aAr_i) (5.7) 

Q<i<N-l 

Since when a particular a^ vanishes the variable Xi only appears with non negative power in 
V{xq, ...xn-1', Go, ■■■, «i = 0, ...ttN-i) we can let x^ = in the computation of the constant 
term, i.e. 

G'Ar(ao, ■ ■ ■ , Oi = 0, ■ ■ ■ , aN-i) = Gn-i{clo, ■ ■ -, di, ■ ■ ■ , aN-i) (5.8) 

Furthermore 

G'Ar(0,---,0) = l (5.9) 



Finally Gn is obviously symmetric in the a^'s. This property together with ( |5.7[ ) ( |5.8D and 



uniquely determines Gn recursively and it is clear that the r.h.s. in (|5.5| ) satisfies 
the same conditions. This suffices to prove the required equality. 

Let us apply Dyson's theorem in the following form. The integral 

27r 2n 



is equal for p = 2s + 1 odd to 



^' E \slt,.-f (5-11) 

0<Zo<---<Ziv_i<p(Af-l) 
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according to its expansion ( 2.11 ). On the other hand it is also equal to 



2n 



2n e^(^o+---+^^-i)(^- 



^(-1)' 



duQ duN—i 

~2^ 



27r 



n 



0<A;<JV-1 



f,i{N-l)Bk 



(5.12) 



{Np )\ 



\\N 



where the last equality is obtained from ( |5.5D by setting all a^ = p. Therefore for p odd 
we conclude that we have the sum rule, 
Property 6 

/ . lfi'Zn,---J«_i I = ^^TTTZnTv (5.13) 



'0,'",'iV-l I 



0<Zo<---<Zjv-i<p(Af-l) 



N\{j>\)^ 



where both sides are integers (i.e. A^! divides t-Mt). For example when s=l we have 



{Nv)\ 



{p\y 

El 

0<Zo<---<iiv-i<3(Ar-l) 



(SAT)! 



,(^) 



In a similar manner we can compute the coefficient g^,^_Y. gifq_x\^\ (s^x^iN-V) 
most compact term noting that the corresponding character is o"^ ~- 

hence 

TV' o^''^ - 

dB^ d6N-i A(e^^°, ■ ■ ■ , e^^^-i)2^+iA(e-^^°, ■ ■ ■ , e 
~2^ 



(llo<fc 



-Wn-i' 



(5.14) 



of the 



<N-l-^k I 



27r 



no 



<A:<Ar-l 



,iefe(Af-i) 



(5.15) 



, N(N-x)^ r d9o d0N-i -TT 

fc=0 



^-lrp,(^.,,).s+l 



27r 
This gives 

fi's(Ar-l),--,(s+l)(iV-l) 

In particular when s = 1, 

gN-l,N,---,2N-2 = 



2-K 



,i{N-i)ek 



-1 



NjN-Ds [(s+l)A^]! 



iV!(s + l)!^ 



Ar(JV-l) 

-1)^^-^(2A^-1)!! 



(5.16) 



(5.17) 



From equation ( |2.17| ) and ( ^.13|) it follows that, 

^^^^- inf 1,\...1^_,\<Zn{p)<^-^ sup /o!- 



(p!)^ Yadr 



(pl 



IW 



■In-i\ 



(5.18) 



Yadm 
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The sup and inf are taken over admissible tableaux. We suspect that these correspond 
respectively to the most extended and most compact terms, in other words that 

m n ^s^.w^n^ (,19) 

^^'' k=s{N-l) ^' 

Referring to equation (A. 18) and our expectation ( |2.32| ) we see that these bounds are 
unfortunately quite miserable. We add the following remark. Consider for simplicity in 
Z]^ (3) the ratio u^ of the contribution of the most extended to the most compact term, 

nS"-i^!((2iV-l)!!P 

Then, 

UN+i _ , iV!(3iV)! ^3_ 27 ;v ., .^^ 

~^ ~ ^2iV + l)!2 ~^^- (2iV + l)2^16^ ^^-^^^ 

Therefore for N —>■ oo, un goes to infinity which suggests that in the large A^ limit, 
Zjv(3) is dominated by the contributions from extended states in agreement with "physical 
intuition" . 



6. Linear system for the coefficients 

We now come to linear equations satisfied by the expansion coefficients g, where for 
simplicity we will set s = 1 although all that will be said admits a generalization to arbitrary 
s. We start with the following obvious observation. Let Q{xo, ■ ■ ■ , xn-i) be an arbitrary 
symmetric polynomial in xq, ■ ■ ■ , xn-i, homogeneous in these variables of degree A^(A^ — 1) 
and vanishing when xq = xi. Then Q{xo, ■ ■ ■ , xn-i) is proportional to A(xo, ■ ■ ■ , xn-i)^- 
Indeed it admits the factor xq — xi, hence from symmetry is divisible by A. But ^ is 
antisymmetric, thus is divisible by A. From the comparison of degrees, -^ is a constant. 
Similarly an antisymmetric polynomial in xq, ■ ■ ■ , xat-i, of degree ^N{N — 1) vanishes 
obviously at xq = xi. If its derivatives also vanish at xq = xi then it is proportional to 
A(a;o, ■ ■ ■ , xn-i)^. Let us apply this to the expansion ( |2.1C1| ) - ( p.ll| ). By requiring that 
the sum 

J2 gi,,..,i^_Ax'"---x'--^\ (6.1) 

0<io<---<'jV-l<3(^-l) 
S!j=3JV(JV-l)/2 
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obviously antisymmetric in xq, ■ ■ ■ , xn-i has vanishing derivative at xq = xi, we will obtain 
a quantity proportional to A^. The condition reads (by multiplying it by xq assumed 
different from zero), 



o<;o<-'-<'jv-i<3(Jv-i) 

Sii=3iV{iV-l)/2 



yio,'",lN-l 



Jo 



IqXq I-iXq 



X 



X 



Jur 



X 



h 






I 



N-1Xq 



X 







X 



Zjv-i 
N-1 



Let us expand this in powers of xq with coefficients which are determinants in X2, 
We set equal to zero the coefficient of 



(6.2) 



,xn-i- 



X 



la+lb \^l 





X 



. X^o. . . . x^b . . ■x\{x2, ■ ■ ■ , XN-i) 



extending gi^^...^ij^_^ as an antisymmetric tensor in its indices to get 
Property 7 

Given non negative integers < /q < ■ ■ ■ < In-^ < 3(A^ — 1) with sum equal to 

— L ~ ' — L, L > 0, the coefficients g. extended as antisymmetric tensors are uniquely 



defined by the conditions 
Y^ {bM))gi^^...^i^_^^i>^i=0 

(6.36) ^o,3,-,3(iv-i) = 1- 



l)>!>i'>0 



Examples: 

N = 2 
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l>l'>0 
l + l'=3 



{I -l')9i',i = ^90,3 + 91,2 = (6.4) 



(6.5) 



i-e. 90,3 = 1 and gi^2 = -3. 

(a) ^90,3,6 + 90,4,5 = 

(b) 4^1,2,6 + 2^1,3,5 = 

(c) ^92,1,6 + 92,3,4 
= -%1,2,6 +92,3,4 = 

(c?) GfiTg^o.e + 45r3^i^5 + 2gs,2,4 

= -6fi'0,3,6 - 4firi^3^5 - 25f2,3,4 = 
(e) 55f4,o,5 + 94,2,3 

= -%0,4,5 + 92,3,4 = 
(/) ^95,0,4 + '2.95,1,3 

= 46-0,4,5 + 2firi^3^5 = 
{9) "^96,0,3 + 96,1,2 

= 3(70,3,6 + 6*1, 2, 6 = 

Conditions (a) and (g) yield from the normalization 

6*0, 3, 6 = 1, 6*0,4,5 = 6*1, 2, 6 = —3 

Conditions (b) and (c) give the remaining coefficients 

fi'1,3,5 = 6 , 92,3,4 = -15 

One checks that (d), (e) and (f) are identically fulfilled. This example illustrates the fact 
that we have an overdetermined set of equations. To make the algorithm really useful 
would require to extract a subset of equations of co-rank equal to 1. This means that if 
An = Aj^ is the number of admissible Young tableaux to find a subset of An equations 
of rank An — 1. Should this be possible we could solve for the (7's in terms of minors of 
the corresponding matrix up to normalization. 
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A weak consequence follows from the fact that A(xo, ■ ■ ■ ,XAr_i)^ vanishes when all 
Xi = 1 hence, since c/iy (1, ■■■,!) = dirriY- 

J2 9YdimY= Yl 9io,ii,-,iN-iY[lj:rf]=^ (^'^^ 

Y adm 0<Io<---<In-i<3{N-1) i>j ^ ^ 

where dimy is the dimension of the corresponding representation of the linear group. 
When N = 2 this coincides with our previous condition, while for A^ = 3 one finds 

= 27-10x3-10x3 + 8x6-15 = 
This may be interpreted as saying that A(xo, ..., xn-i)'^ is the difference of two characters 
of GLn operating in two vector spaces of equal dimension 

y qy dimy 

Y adm,gY>0 

7. Specializations 

In this section we study specializations of the discriminant leading to sum rules on 
the coefficients gy ■ For simplicity we stick again to the case s=l. 
As we saw when Xk = e^^ , 

a9/. 2i7r 2i7r(iV-l) ^ (iV-l)(iV-2) .^ ,„,, 

A2(l,e~,---,e — ^ — ) = (-!) ^ N^ (7.1) 



For any integer / let < / > denote its representative mod N in the range 0, 1, • ■ • , A— 1. If the 
sequence < /q < ■ • ■ < In-i < 3(A — 1) is such that in the sequence < /q >,■■■, < In-i > 
two numbers coincide, the corresponding character c/izq,...,;^_j(1, e~iv-^ . . .) vanishes. If all 



< 


lo>,-- 


■, < 


In- 


_i > are 


distinct let e(/o, 

f<lo> 
V 




', In- 
< 


-l) 

In- 

A- 


denote the 

-1 ; 


sign 


L of the permutation 

(7.2) 


in 


which 


case 


we 


have 
























chi„,- 


, (1 e N ... 


, e 


N 


^) 


= e(/o,---, 


In- 


i). (7.3) 



Thus, 

(Af-l)(iV-2) 



iV^(-l) "'^ = J2 9io,-,i.-.<lo.---.lN-i) (7.4) 



y adm 
{<lo>---<l'N~l>}p^'^'>^- °f {0---N-1} 



46 



The reader can check this sum rule using the hstings in appendix D. 

More generaUy we can speciahze the x^'s to be the successive powers of a single variable 
g, the so called principal specialization, 

x^ = q' Q<i<N -1 (7.5) 

Thus 

0<i<j<Ar-l (7.6) 

= (-i,^^,Er-"(— ).[iy2],...[iv-i], 

where the g-factorial is defined as 



[n], = {l-q){l-q')---{l-q-); Inn -i^ = n! (7.7) 



and 

7V-1 



A^(A^-l)(A^-2) 



E^(^-l-^)= ^ 6 ^ ^ (7.8) 

n=0 

Note that the symbol [oo\q makes sense for |g| < 1 being equal to q~^''^^rj{q) where rj is 
Dedekind's function, and that Jacobi's identity reads 



n>0 j>0 



Under the same specialization. 



0<a,b<N-l 

N(N-l) 



0<i<j<N-l 
0<i<j<iV-l 



(7.10) 



It follows from (Q and (^) that 

o<io<--'<'iv-i<(2a+i)(iv-i) 0<i<j<Ar-l 

(7.11) 
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Set 

Hk = lN-i-k-i2s + l){N-l-k) 0<k<N-l (7.12) 

which represents, counted from the "top" the departures of the Vs from their values in the 
most extended case. The exponent in the expUcit g-factor on the r.h.s. of ( [7.10| ) reads 

N-l N-l N-l 

E{lo, ■■■, In-i) =J2^h- {2s + l)k]{N - 1 - A;) = ^ fc^Ufc = ^ nfc > (7.13) 

fc=l k=l 



in terms of the notation introduced in section 4. Its extreme values are 



E{0, 2s + 1, 2{2s + 1), • • ■ , (AT - l)(2s + 1)) = 



(7.14) 



fe=l 



Since the power of (1 — q') occurring on the l.h.s. is higher than the one in any term on 
the r.h.s. we recover ( |5.(j| ) and its generalization to all s by letting q' ^ 1 on both sides 
of (|7.11| ). Of course one can equate further terms to zero a,s q —>■ 1 (in fact as many as 
sN{N — 1) ) to require a similar leading behavior on both sides, but this would be less 
effective than the equations discussed in section 6. 
Finally we can rewrite (|7.11|) as 



([i],---[iv-i].r= E ^i^^^?^-^ "^^"^^ n 



Y admiss bEY 



1 - qHb) 



(7.15) 



where following Macdonald b runs over the boxes of a tableau, h{b) is the hook length 
at 6, i.e. one plus the number of boxes at the right or below 6, and c(6) is its "label" or 
"content", i.e. the sum i + j where i is the line index starting with zero from the A^-th 
bottom line and j its column index starting at 1 at the left as shown in the illustration for 
A^ = 4 where the shaded box has h{b) = 7 and c{b) = 6. Equivalently the label of the left 
uppermost box is A^ and labels increase (decrease) by a unit when moving one step to the 
right (down). 

Fig. 2: Hook length and label of a box 
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n s^ = ^^"^^(^^A^) (7-18) 



The expressions occurring on the right hand side of equation ( [7.15|) generahze the 
g-binomial coefficients. For a tableau with a unique column Y^ of height k < N we have, 

n^[i_qmj (i_g)(i_g2)...(i_g.) [kUN-k], ^-^ 

while for a single row of length k, call it Y^, 

TJ f ^-^'^'^ \ (l-9^)(l-g^+^)---(l-9^+^-^) [N + k-1], 
beVfc 
and as g ^ 1, they reduce to 

c{b) 
beY 

which gives one of the fastest means to compute these dimensions. Even though equations 
( [7.11| ) or ( |7.15| ) depend on a single variable q it is still a tremendous task to extract from 
these equalities the coefficients g for instance by expanding at small q. 



8. The number of admissible tableaux 

The number of admissible tableaux A)^ is according to property 3 of section 4. the 
number of A^ + 1-uples of non-negative integers (no = 0, ni > 0, ...,nN-i > 0, tt-at = 0) 
such that 

2s-2nk+nk-i+nk+i>0 /c = 1, 2, .., A^ - 1. (8.1) 



Fig. 3: The polytope 114 and its 16 integer points. The black dots belong 
to vertices or edges, the half-filled to faces, and the empty one is strictly 
inside the polytope. 

(s) 

This counts the number of points with integral coordinates in the closed polytope nj^ 
defined by the equations (|8.1|). For illustration, we represent the 3-dimensional polytope 
7r4 and its 16 points with integer coordinates on Fig. 3. We first concentrate on the number 

An = A''^\ for s = 1. 

Let /Ar(a;o, ...,xn) be the Laurent series 

22 2 

Jn[Xq,...,Xn) 
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and expand its denominators as infinite power series to find 

N-l 



/-= E n 



2 ( ^j-l^j+1 
^J 1 ^2- 



ni,...,njv_i>0 j = l \ J 

N-l 



E^ni nM-\ TT 2-2nfc+nfe_i+nfc+i 
"1, ■.,"jv_i>o A;=l 

tlierefore the desired number A^ is equal to tlie number of monomials appearing in [/Ar] + , 
the polynomial part of /at. 



Examples 

1) N=2: 
We have 



X0X2 • 



and 



[f2{xo, Xi, X2)]+ =xl+ XqX2, 



thus A2 = 2. 










2) N=3: 










We have 


/3(xo,Xi,a:2,X3) = 




X^X2 






(1- 


-B^){1- 


XIX2, \ 



and 



\f-i{xQ, Xx,X2-, 3^3)]+ = 2^13^2 + ^03^2 + ^S^^l + 3^03:1X20:3 + x\x\, 



which implies A3 = 5. 

3) N=4: 
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T^ _|_ — iL-i JunJurt '\ tZ/ Q tZ/ n X Q ~i~ X 1 1// Q ~i~ X -| Jb'^Jb/^ n~ XqX]^X2Xq ~i~ X q X q 

i4_i_3 _i_ 2 _i_22_i_22 

"T" XqXoX^4 ~r =^1 X2t^3'^4 "r XqX]^XoX3'^4 "t XqXi XqX^4 "t XrjX2'^Q'^4 

i42| 2 2|222|2 2,33 

hence A4 = 16. 

We now derive an integral formula for A^r. Let / be a Laurent series of the form 

m 

fc = — CXD 

and /+ = XlfcLo '^fc'^'^- -^y Cauchy's theorem, we can write 

i|x|>i2z7rx-l ^ 

This gives an integral formula for A2, by replacing / by /2., and taking xq = X2 = ^- More 
generally, we have a similar formula for 

A^ = [/Ar] + (1,1,...,1), 

which reads 

, (iXi (iXAT-i /Ar(l,Xi,..,XAr_i,l) 

An — 



for any p > 1. The choice of contours ensures that 

I Xk-lXk+l 



2 



<1 and \xk\>l k = 1,...,N -1. 



The number of states Aj^ corresponding to the expansion of A^* is obtained by usin^ 
as generating function 

„{s), N _ -t-1 J-2 ■■■-t-AT-l 

J^ [Xo, Xi, .., XtvJ - ^_i Xfe„iXfc + i ^ ' 

11a;=1 U ^2 J 



and 



O^- /-I IT' - -^ 



dxi dxiv-i /^^^(l,Xi,..,XAr_i,l) 



|xi|=p»(~-o/2. 2z7r 2z7r 11^=1 (a^j - 1) 
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for arbitrary p > 1. The computation of A]^ is consequently reduced to a straightforward 
but tedious extraction of residues. The result for A]^ is a polynomial of degree (A^ — 1) 
in s with rational coefficients. This is a general property as discussed below. 

We list below the first few polynomials A^^. 

a['^ = 1 

A^^^ =s + l 

aI'^ = 2s^ + 2s + 1 
^(s) 16 o 13 2 19 

■*" = ¥' +y'+ir'+i 

a(s) 50 4 , o 40 o 
^3 3 

As) _ 

,(5) 118098 o 279936 7 367144 ^ ^ . 75724 , 5911 . 

^ 35 35 45 45 15 

19927 2 163 

H s H s + 1 

315 21 

4(-) _ 
^10 ~ 



^1 



3 2 6 

50 4 . s 40 2 . 

288 5 ^^ 4 385 3 ^^ 2 157 
: ^s^ + 96s^ + —^s^ + 23s^ + —-S + 1 
5 6 30 

9604 fi 6076 . 2858 4 _, o 1531 2 267 

: s^ H s^ H s^ + 134s^ H s^ H s + 1. 

45 15 9 45 45 

262144 7 26624 . 71992 . 6371 4 43657 . 5783 2 971 

315 15 45 8 180 120 140 
118098 = 279936 . 367144 « .„_„ . 75724 . 5911 . 



s H s + i 

315 21 

8000000 n 2320000 o 39666608 7 1236328 . 12340889 

567 63 945 45 1080 
228025 4 13641133 o 204457 2 10957 

72 22680 2520 1260 
js) _ 857435524 ^q 162983612 9 205427098 § 50279276 7 

^^ 14175 945 945 315 

51136132 ^ 365944 . 61904509 4 1628849 

675 15 11340 1890 
622849 2 323 

6300 35 



A^: = 



-s' 



(8.2) 



These quantities take a neat form when expressed in terms of the polynomials 

s{s-l){s-2)...{s-k + l) 



CTkis) 



k\ 
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instead of powers s^. We have 

a['^ = 1 

4'^ = (71 + 1 

4'^ = 4(72 +4(71 + 1 

A['^ = 32(73 + 45(72 + 15(71 + 1 

A^""^ = 400(74 + 744(73 + 404(72 + 58ai + 1 

A|;^^ = 6912(75 + 16128(74 + 12481(73 + 3503a2 + 246ai + 1 

A^^^ = 153664(76 + 432768(75 + 437776(74 + 188244(73 + 30702(72 + lllOai + 1 

A^""^ = 4194304(77 + 13860864(76 + 17367872(75 + 10162473(74 + 2731521(73 

+ 275599(72 + 5301ai + 1 
aI'^ = 136048896a8 + 516481920a7 + 773038912(76 + 574771360(75 + 218829232a4 

+ 39174798(73 + 2537596a2 + 26375(7i + 1 
aI'J = 5120000000(79 + 21964800000(78 + 38261688576(77 + 34587246976(76 

+ 17164897361(75 + 4532740981(74 + 561999521(73 + 23932596a2 + 135669(7i + 1 
A['^ = 219503494144(710 + 1050351430656(79 + 2088303502080a8 + 2225685070528(77 

+ 1366844046336a6 + 482748121856(75 + 92013182474a4 

+ 8109500738a3 + 230617119(72 + 716541ai + 1. 



(8.3) 

These triangular expressions for A)^ in terms of the (Tk{s) are easily inverted as 
follows. If 



N-l 
k=0 

then one has 



a'N,k = 

s=0 



Y,<Qi-r>' 



due to the orthogonality relation: 



EC !)M)'-"-*.... 



fc=Tn 



This shows why the coefficients in (8.3) are integers. We have no clue as to why they 
turn out to be positive, implying inequalities among A]^^ for fixed N. 
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The problem of counting the number of integral points in an integral polytope (i.e. a 
polytope with all its vertices at integer points) is a classical (and difficult) one. We base 
the following on the work of Ehrhart [0, who, starting from an integral (i-dimensional 
polytope n with vertices in 2^ considered the number of lattice (integral) points in the 
dilated polytope sll, s integer, 

Eu{s) =card(Z'^nsn). 

He showed that this number is a polynomial of degree din s, hence deserves to be called the 
Ehrhart polynomial of H. Our number A)^ of admissible tableaux is precisely the Ehrhart 
polynomial of the (A^ — l)-dimensional polytope 11]^, defined by the equations (P?T|), with 



•-N ' 

Hi G lR-(-. In the following, we derive a simple expression for the Ehrhart polynomial of 

a specific class of polytopes, which we conjecture is valid for 11]^ • We will make use of 

the Ehrhart reciprocity theorem, which reads as follows. Let -E'n(s) denote the number 

of lattice points which lie strictly inside the dilated polytope sll, i.e. the number Eu{s) 

minus the number of lattice points on the boundary d{sll). Then one has the reciprocity 

relation 

Theorem (Ehrhart): 

Eu{s) = i-l)''Eui-s). (8.4) 

Given an integral polytope U in d dimensions, we wish to perform a decomposition 
into "elementary cells" , namely non-flat polyhedral integral simplices with d+ 1 vertices, 
which do not contain any other lattice point (the intersection of an elementary cell with 
Z'^ is reduced to its ci + 1 vertices). We have the following 

Proposition 1. Such a decomposition exists but is in general not unique. 
To prove the existence, we rely on a classical theorem which guarantees the existence of a 
decomposition in any integral polytope. Let us consider a simplex in this decomposition. If 
it is not an elementary cell, it contains at least one lattice point P. If we draw all the edges 
linking this point P to the (A^ + 1) vertices of the simplex, this defines a decomposition of 
the simplex into smaller simplices. 

Fig. 4: The three possible new decompositions of a tetrahedron, according to 
the position of the point P on an edge (2 tetrahedra), on a face (3 tetrahedra) 
or strictly inside (4 tetrahedra). 
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Retaining only the non-flat ones, we end up with {d + 1) simplices if the point P is 
strictly inside the initial simplex, d simplices if it lies strictly inside a {d — 1) dimensional 
face,..., 2 simplices if it lies strictly inside a 1 dimensional edge of the initial simplex 
(we illustrate the case d = 3 on Fig. 4). We can iterate this procedure with the new 
decomposition, until the decomposition is into elementary cells. 

For any such elementary cell decomposition of H, the total numbers of vertices (0- 
faces), edges (1-faces),..., /c-faces,..., cells (li-faces) are not in general invariant. However, 
the situation gets much better in the following particular case 

Proposition 2. If the polytope 11 is such that one of its elementary cell decomposi- 
tions is only made of "basic cells" , i.e. elementary cells with volume equal to the minimal 
value ^, then the numbers of j-faces of the decomposition, which we denote by Fq, -Fi,..., 
Fd, are arithmetical invariants of the polytope H. Such a polytope will be called basic in 
the following. 

To prove the above, we will derive a formula for the Ehrhart polynomial of H, involving 
only the Fj, j = 0, ...,(i, which will consequently appear as invariants of H. The issue of 
whether an elementary cell is basic or not becomes relevant in 3 or more dimensions. For 
d = 2, an elementary cell is a triangle with vertices in Z^, which contains no other lattice 
point. Any two of its edges define two basis vectors of Z^, giving rise to a parallelogram of 
area 1, therefore the triangle, half of the parallelogram, has area |. So in 2 dimensions, all 
elementary cells are basic. This is no longer the case in 3 dimensions. Take for instance the 
tetrahedron with vertices (0,0,0), (1,0,0), (0,1,0), (l,p, q'), p and q co-prime, it defines 
an elementary cell, with volume |, in general not basic. Nevertheless one checks that the 
closed tetrahedron intersects the lattice Z"^ only at its vertices! However, in the particular 
case of H4 ( of Fig. 3 ), it is possible to show that it decomposes into 32 elementary cells 
(tetrahedra) which turn out to be all basic (volume |). We are thus led to the conjecture 
that the polytopes Hjy are basic, for all A^. 

Before treating the general case of a basic polytope, let us first concentrate on the 
simplest example of the basic simplex E^, with vertices at (ni, ..., n^), where either all n^ = 
0, or all Hi = except one, which is 1. The simplex itself is its own basic decomposition, 
and we compute very easily the numbers F- of j-faces 
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On the other hand, the computation of the Ehrhart polynomial Ey,^{s) 
is very simple, we find 



card(sEd n Z° 



E.jA = E E 



1. 



A;=0 "l,"2.-'-."d>0 



The number of partitions of the integer k into d arbitrary integers is just ( ^ ) , and 



fc=0 ^ ^ 



{s + l){s + 2){s + 3)...{s + d) 
d\ 



which can be recast into 



->'^ge:D( 



d+l\ fs-1 

j 



In this particular case, we get a simple expression for the Ehrhart polynomial in terms of 



theF, 



(d) 






s-1 

j 



(8.5) 



We are now ready to generalize this to any basic poly tope 11 . The result reads 
Proposition 3. Let 11 be a basic (i-dimensional integral polytope, and Fj, j = 

0, 1, ..., (i denote the numbers of j-faces in a given basic decomposition. Then the Ehrhart 

polynomial of 11 reads 




(8.6) 

Proposition 2 above appears as a simple corollary: due to the invariant definition 
of Eyi{s), the numbers Fj are also invariant, i.e. do not depend on the particular basic 
decomposition performed. Note also that the intuitive fact that Eii{s = 0) = 1 for any 
polytope n, because the zero-dilated polytope is reduced to a point, translates into the 
non-trivial Euler's relation linking the numbers Fj to the zero genus of the (i-dimensional 

decomposed space 

d 

J2(-^yFj = ^n(O) = 1. 
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Let us proceed to the proof of proposition 3. Equation ( |8.6| ) was just derived in the case of 
the particular basic simplex E^^, and extends trivially to any elementary basic cell. For a 
given polytope 11, we first perform an elementary cell decomposition as in sect.l. We want 
to count the number Eyi{s) of lattice points in the dilated polytope sll. This number is 
equal to a sum of (d + 1) terms. 

(0) the number of dilated vertices, Fq. 

(1) the number of dilated (1 dimensional) edges of the initial elementary basic cell 
decomposition, Fi, multiplied by the number of lattice points which lie strictly inside each 
dilated edge. This last number is (s — 1), the same for all the dilated edges. 

(j) the number of dilated j-dimensional faces of the initial elementary basic cell de- 
composition, Fj , multiplied by the number of lattice points which lie strictly inside each 
dilated j'-face. By Ehrhart reciprocity ( |8.4| ), the number of lattice points inside any dilated 

j-face is 

^^ {-s + l){-s + 2)...{-s + j) 



-lVi?^,(-.) = (-ll 



Jl 



(s-l)is-2)...{s~j) 



j! 



s-1 

J 



[d) the number of dilated d-dimensional basic cells of the elementary basic decom- 
position, Fd, multiplied by the number of lattice points strictly inside each cell, namely 

(s-l)(s-2)...(s-N) 

Summing all these contributions, we get the desired result ( |8.6| ) 

j=o \ -^ / 

When applied to our case of (conjecturally) basic polytopes HJy , the above formula 
( p.6|) suggests to rewrite the Ehrhart polynomials A)^ in the basis 

{s-l){s-2)...is-k) fs-l\ 

instead of crfc(s) as in (8.3). We get again positive integers, which reinforce our hope that 
our conjecture might be true. We obtain the following table 
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a['^ = 1 

Ai'^ =1,1 + 2 

A^'^ = 4z/2 + 8z/i + 5 

A["^ = 322^3 + 77z/2 + 60z/i + 16 

aI"^ = 400/^4 + 1U4U3 + 1U8U2 + 462/^1 + 59 

A'^q^ = 6912Z/5 + 23040Z/4 + 28609z/3 + 15984z/2 + 3749z/i + 247 

A'f^ = 153664Z/6 + 586432Z/5 + 870544z/4 + 626020z/3 + 218946z/2 + 31812z/i + 1111 

A'^"^ = 4194304/^7 + 18055168/^6 + 31228736/^5 + 27530345/^4 + 12893994/^3 

+ 3007120//2 + 280900//1 + 5302 
A^^^ = 136048896/^8 + 652530816//7 + 1289520832/^6 + 1347810272//5 + 793600592/^4 

+ 258004030/^3 + 41712394/^2 + 2563971/^1 + 26376 
A["J = 5120000000/^9 + 27084800000//8 + 60226488576/^7 + 72848935552/^6 

+ 51752144337/^5 + 21697638342/^4 + 5094740502//3 + 585932117/^2 + 24068265//i + 135670 
A[\^ = 219503494144/^10 + 1269854924800/^9 + 3138654932736/^8 + 4313988572608/^7 

+ 3592529116864/^6 + 1849592168192/^5 + 574761304330/^4 

+ 100122683212/^3 + 8340117857/^2 + 231333660/^1 + 716542. 

(8.7) 

For A^ = 3, we read that the 2-diniensional polytope Ilg has a basic decomposition 
into 4 basic triangles, with 8 edges and 5 vertices. For A^ = 4, we have 32 basic tetrahedra, 
with 77 triangular faces, 60 edges, and 16 vertices, etc... 

In particular, if our conjecture is true, the volume of the polytope 11]^ should read 



Vn 






(A^-l)!' 



where -Fat-i is the number of basic cells of any basic decomposition of 11]^. In the remain- 
der of this section, we will derive explicitly the first few leading coefficients of A]^ when 
s is large. The leading coefficient is intuitively the volume V^ = ^ (n-i)\ — (see appendix 
E for the calculation) of the polytope IIJy ■ A more surprising feature is that more can be 
said on the subleading coefficients. 

We apply a particular case of the combinatorial Riemann-Roch theorem of Kantor 
and Khovanski |jT^ . The latter relates the number of points with integer coordinates in 
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a given integral polytope to the volume of the polytope deformed by slightly shifting its 
faces, parallelly to themselves. To simplify notations it proves more convenient to work in 
A^ dimensions, i.e. to consider the polytope Hj^.].! and its Ehrhart polynomial A^j^^. 

Let yjY'_|]^(/ii, ..., /lyv; ci, ■■■, CAf) denote the volume of the deformed polytope defined 
by the inequalities 

-2efc < 2nk < rik-i + rik+i + 2s + hk k = 1, ..., A^, 

with as usual no = un+i = 00, and T{x) the generating function of Bernoulli numbers 

T(x) = = 1 + -X H x^ x^ H x^ + ... 

^ ^ 1-e-^ 2 12 720 30240 

Consider the polynomial Bn{s) of degree A^ in s 



Bn{s) = 






'-fc=i ^^^ ^^'^ 



hi=£i=0 k=0 






On the other hand, the desired number of points, A^_^^, is also a polynomial in s 



N 



^i^E"!"'""- (8.9) 



fc=0 

The theorem relates the coefficients of both polynomials, we have 

(N) _ , (TV) _ .. 
(N) AN) 

(N) _ , (TV) 



LIO) 



^N-2 — "JV-2 



bN-2+ Yl f^N-2{FN-2)rk{F2), 



where Vn+i is the volume of the polytope n^:!|_j^, the last sum extends over the (A^ — 2)- 
dimensional faces Fn_2 of the polytope HJ^^j^, jjn-i is the so-called (A^ — l)-dimensional 
relative measure of n]^:, j^ (the sum of the (A^ — l)-dimensional measures of the (A^ — 1)- 
dimensional faces of 11^^!^^, such that the unit cell of the restriction of 2^ to the given 
face be taken as unity), fXN-2{FN-2) the relative (A^ — 2)-dimensional measure of the face 
Fjv_2 defined analogously, and T2 a certain invariant of the two-dimensional polar cone 



This definition is only valid for A'' > 1. When N = 1, one has to take the deformation 
-ei < rii < s + hi, in which hi is multiplied by 2 in the definition. 
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^2 to Fn-2 (the cone generated by the normals to the two (A^ — 1) -dimensional faces of 
which Fn-2 is the intersection). 

The crucial hypothesis which enables to apply the above mentioned theorem is that 
the polytopes we are considering have so-called 2-primitive fans. The fan of a polytope is 
the set of cones generated by subsets of the (integer valued) normal vectors to its (A^ — 1)- 
dimensional faces: it is a kind of dual envelope of the polytope. The /c-primitivity means 
that for all rn-dimensional cones of the fan, m < k, the generators of the cone, i.e. a set 
of m (integer valued) normals to the (A^ — 1) -dimensional faces of the polytope, can be 
completed into a 2-basis of the lattice Z^. In our particular case, the normals to the 2N 
{N — 1) -dimensional faces read 



■th 



(0, ..,0,-1,0, ...,0) - 1 in i'" position i = 1 
(0, .., -1, 2, -1, 0, .., 0) 2 in z*'' position i = 1, .., A^ 

It is clear that for instance the two normals 



A^ 



(2, 



-1,0, 



0) and (0, -1,0, ..,0) 



cannot be completed into a Z-basis of Z . Otherwise one could find (A^ — 2) vectors of 



Z-^ fi, f2, ••, ^Af-2, with integer valued coordinates v. 
(volume of the unit cell associated to the basis) 



(J) 



G Z, such that the determinant 



2 

-1 -1 





,(3) 







vl 



(N) 



^N-2 
^N-2 



^N-2 



would be equal to one. But after subtracting the second column from the first one, we 
can factor out 2, hence the determinant cannot be equal to 1. Therefore our polytopes are 
2-primitive, as any single normal can be completed into a Z-basis of Z-^. 

The above proves rigorously the intuitively obvious fact that the volume Vat+i of the 
polytope n^^j^ appears as leading s^ term in A'^'_^^, and gives us compact formulae for 
the next two subleading coefficients. 

We summarize the conjectures in the following proposition, and leave the detailed 
arguments to appendix E. 
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Proposition 4. For large s the polynomial A^,^ behaves like 

An+i{s) = a^^h^ + a^^\s^-' + a'^^l.s^-' + 0{s^-% (8.11) 

where 

a^ -Vn+1- ^j 

N 1 ^ 

k=0 



(N) _.{N) 1.. 

«7V-2 -»Af-2 + 4'^^'-l 



2Af-2 

■ [(A^ + 1)^-^(3A^^ + 17A^^ + 72N^ + U4N + 266) + 6(A^ - 1)^"^] , 



(A^-2)!4! 

(8.12) 

Our conjecture that n]y;_|^j^ is basic would imply that the various numbers F^r, F^^i 

and Fn-2 of respectively basic cells, (A^ — l)-faces and (A^ — 2)-faces in the basic cell 

decomposition of n]Y_j_i, are given by 

Fn-1 = 2^-^{N + 1)^-^(2A^^ + 7N^ + 9N+ 10) 

ON-A (8.13) 

Fn_2 = \(N + 1)^-^(12N^ + 64N^ + U3N^ 

6 

+ 237A^^ + 264A^2 + 212A^ + 258) + 6(A^ - 1)^"^] . 
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Appendix A. Asymptotic evaluation of lnZ]\[{l) using Euler- Mac Laurin's for- 
mula 



As indicated in eqn.(p.7|) we have 

AT 

Zn{1) = 1[3'. (A.l) 
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This section is devoted to a proof that asymptotically as A^ — i> oo 

lnZN{l) = InN - -N'^ + NlnN + N{-ln2TX - 1) + —InN 

2 4 2 12 CA 2) 

l-7 + 5/n27r C'(2) 1 1 1 ^, 1 x 

12 27r2 12A^ 720A^2 360A^3 ^ v^4^ 

Here 7 is Euler's constant, C,{t) = X]n>i ^~* ^^ Riemann's zeta-function, and the numerical 
constant is 

l^2±l!^ + m = .75351738... (A.3) 

12 27r2 ^ ^ 



Such an expression vindicates the estimate (|2.32| ) with 



-2/(l) = ^/n27r-l. (A.4) 



Note that beyond the term of order A^ we find corrections of order InN (and not v -^ as 
crudely expected from the thermodynamic reasoning). We guess that this is not peculiar 
to the case p = 1. 

We rederive Euler-Mac Laurin's formula for our purpose []TT|. Let f{z) be an analytic 
function for Re{z) > such that in this half plane \f{x + iy)\ = o{e^'^y). We have in mind 
a function such as f{z) =2", a real positive, and we take the principal determination. 

Fig. 5: The contour of integration for Euler-Mac Laurin's formula. 
Remarking that in the vicinity of any integer /c, cotanyr^ behaves as 

l/7i{z-k) 

we have for a sum /^r. 

In = /(I) + /(2) + ... + f{N) = 1 / dzf{z)cots.n7rz, (A.5) 

2« Jc 

where the contour C is indicated on figure 5, 9 being an arbitrary real parameter between 
and 1 and F a half circle of arbitrary small radius p. Call C_(_ and C_ the upper and lower 
parts of the contour C, excluding the arc F (see fig. 5 ). Noting that for any fixed x 

lim cotan7r(x + zy) = - sgn{y), 

\y\^oo I 
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where sgn{y) is the sign of y, we add and subtract the integrals | J^ f{z)dz. Thus 
In = — dzf{z){cota,mrz + i) + — / dzf{z){cota,n'n:z — i) + 

Zt I r> Zt I r' 

1 f If w (^-o) 

+ ^ / dzf{z)cota.n7Tz - - dzf{z) + - dzf{z). 

As p shrinks to zero, the third integral over F tends to ^f{N), while the fourth and fifth 
combine to yield the integral of / along the real axis from 9 to A^. Letting M -^ oo, the 
two horizontal branches in the first and second integral vanish and we are left with 

In = / dxfix) + -f{N) + / dz J}^J 



In the sense of asymptotic series we replace [f{N + iy) — f{N — iy)]/i in the last integral 
by its Taylor series at y = to the effect that 

rN 1 '}(_'\\k roo 2fc+l 

This is indeed an asymptotic expansion provided 

lim f^^+^\N)/f^''\N) = 

N^oo 

\f^^\N + iy)\ < Aj('=)(Ar)e(2--^)l2^l, 

which will be satisfied in our applications. By using a specific function f{z) = e*^ and 
comparing In+i — In — f{N) to its asymptotic expansion, we find for the coefficient of 

/(2fc+l)(JV) 



2(-l)^ f^ ^ y2fc+i 5, 



2fc+2 



(2A; + l)!7o e^'^y-l (2A; + 2V 

where the Bernoulli numbers B^ are defined through 



dy^^tr^ = 7^^^^^ . (A-9) 



?^-E''t (A.10) 

k=o 



and vanish for k odd larger than 1. If /(O) exists, one can let ^ ^ in which case 

C(0)^-i/(0)-/°°f ^"^j:/_'7> . (A.11) 
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On applying eqn.( |A.8| ) to f{z) = z"", a > 0, one obtains 



^ A^«+i A^" _ _^ ^«+i ,,-^ Bo 



?^°-irTT + ^^^"(''>-^^S(^lfT%«<«-i'-("-^*"^"-^'"'(*-i^' 



3=1 k>0 

where 



6'"+^ „ .„. _ . na [^ ^ y" 



CM — = C«(0) = -2 sin— / dy 

a + 1 2 Jo 



e27ry _ I 



-2sin^X] / ^^/^/^e-^'"^^ 



. 7ra ^ r(a + l) (A.13) 

-2sm— 2^ 



2 ^^ (27rm) 

m=l ^ ^ 



a+l 



„ . na T(a + 1) ^, ,. 
= -2sin — -^— — fC(a + l) 

= C(-«), 

since the zeta function satisfies the reciprocity relation 



C(.) = 2r(l - s){2nr-' sin -C(l - s). (A.14) 



The case f{z) = Inz obtained by taking a derivative at a = of the previous estimate 
yields Stirling's formula as well as 9q.Cq,(0)|q.=o = ~C'(0) = \ln2n. Thus 



1 1. - ^ _B2fc+2iV_-2'=-l 

(2A; + l)(2A; + 2) 



faiV!^(iV+-KniV-iV+-fa2. + X: ;^+; (A.15) 



fc>0 



where of course the asymptotic series diverges ! Similarly when f{z) = zlnz we take a 
derivative at a = 1, use r(2) = 1, r'(2) = 1-7, C(2) = ^ and find 



> ?/7i? ~ InN \ InN ^ InN + ^^ 

4^-^ -^2 4 2 12 12 27r2 

^-' (A.16) 



E B2k+2N 
2A;(2A; + l)(2A; + 2) ' 
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Finally from equations (|A.15| ) and ( |A.16| ) 

N 

lnZN{l)=lnl[j\ 

N 

= J2iN + l-j)lnj 



-[y2{Np — q)lnq + >^ k >^ ln{k + rp)] 

g=l 
^ Np 

-ZnJJj!+i?7v(p), 



g=l fc=l r=0 



^ . = 1 



where 



(A.17) 



AT 

= {N + l)lnNl-^jlnj 

2 4 ^ 2 ^ 12 
1 - 7 + 5/n27r C(2) _}___}: ^__ of— 1 

^ 12 ^ 27r2 ^ 12iV 720iV2 360iV3 "^ *^iV4''' 

as stated in equation ( [A.2| ). 

In view of the above remarks, we add here the following estimate 

N 1 ^P I 1 1 

In T\ijp)\ = -lnY\j\ + ^^^Inip^m) - -ln{{Np)\) + OilnN). (A.18) 

The first term on the r.h.s. is -InZ^pi^)- Indeed 

N N-1 p 

In TT(jp)! = 2_. z_.^^ " ^)^^(^ + ''P) 

j=l r=0 k=l 

Np N-1 p 

= -[Npy Inq — y^ y^irp + k — k)ln{k + rp)] 

P g=l r=0 k=l 



(A.19) 



Rn{p) = -'^k'^lnik + rp). (A.20) 

P fc=l r=0 

With 1 < k < p we have 

Af-l AT-l AT 

y ln{rp) < y^ ln{k + rp) < y^ln{rp) (A. 21) 
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thus 

^'-lnip^-\N - 1)!) < Rn{p) < iL^^ln{p^Nl) (A.22) 



P + '^i^r^N-U.r iNn/D r^^ ^ P + ^ i^^^N . 



2 
which yields the required result (|A.18|) . 



Appendix B. Instability of the cuboctahedron configuration 

In the stereographic representation of the sphere on the plane (D, the equilibrium 

dV ^ d\ 
dzi dz 



configurations, stable or unstable, {zj} are solutions of the equations ^ = ^ = 0, Vz, 



that is 

N 

Y. = (Ar-l)-^, i = l,2,...,iV (B.l) 

z_— / r- — r • -4- ^-r- 



j = i 



Z'l Zj 1 ~r Zj^Zi 



expressing that the potential 

V = -Y^ Inz.jZij + {N -1)Y^ ln{l + ZjZj) (B.2) 

is stationary. Since V is invariant under an overall rotation of the configuration, that is 
by a unitary homographic substitution 

_ azj + 7 
%■ ~^ sj - 



-^Zj + a (B.3) 

aci + 77 = 1 

the system ( p.l|) is formally invariant in this substitution. We can easily then show the 
existence of equilibrium configurations on the sphere defined not only by the regular poly- 
hedron but by their composites in diverse polyhedra, so that we have by symmetry 

N-l 

E ^ = (B-4) 

by putting the vertex ^n at the pole {^n = 0). 

Around a solution {zj} of ( |B.1J ) the potential is approximated, up to third order in 
the variations {Szj} by the real quadratic form Q{dz, dz) 

2Q = YJ AijSziSzj + Aij5zi6zj + 2Di6ziSzi (B.5) 
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with 



A 



^j 



dzidzj 



(Z^-Zj) 



+^^. (e 



\ 



^{z,-zir ^^ ^^n^...A2 



{l + ZiZiY 



l^i 



D, 



d^V 



N -1 



(B.6) 



(B.7) 



dzidzi (1 + ZiZi)"^ 

This is a real form in the 2N variables 5xj, dyj {Sz = dx + iSy) with rank not greater than 
2N — 3 since the variations {6zi} associated with an overall rotation of the configuration 
{zi} annihilate Q and depends on 3 real parameters as shown by (|B.3| ). The configuration 
is stable if Q is semi-positive. 

To prove that the configuration 12(24) of the cuboctahedron is unstable we must show 
the existence of a negative eigenvalue of the real 2N x 2N symmetric matrix 



D + ReA -ImA 
-ImA D - ReA 



(B.8) 



but the spectrum of this matrix is not invariant under the rotations. However the signature 
is invariant and this is all we need. To obtain the intrinsic modes of vibration we consider 
the form Q{du, du) 

2Q = y^ Bij6ui5uj + c.c. + 26ui6ui 

with 



B 



1 



^j 



Vd;d, 



-.A 



1-j 



B 



and 



where 



5ui 



iD^dz^ = ViV-1 



j« 



SZi 



(B.9) 
(B.IO) 



4 ^ _ 1 5z5z 
-ouou — -- 



1 + ZiZi 
2 Qji2 



sin^ ed(i)^ + d9' 



(B.ll) 



:b.i2) 



N-l A{l + zzY 

is the invariant metric on the sphere. The secular equation associated with ([B.10|) reads 



BijSuj + (1 — X)dui = 



(B.13) 



and its conjugate. Denoting by p^ the spectrum of the non-negative hermitian matrix BB 
of order A^, we have the 2N eigenvalues 



Xn = I- Pn and 2- Xn = 1+ Pr 
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(B.14) 



and the diagonal form of Q 

N 



Q ^ 5^(1 + pn)den + (1 - Pn)Svl (B.15) 



Tl=l 



With the exception of the three zero modes (pi,2,3 = 1), the configuration is stable if 
p^ < 1, Vn > 3 and it is unstable if max„ p^ > 1. 

The spectrum of BB is invariant under (|B.3| ) since the hessian, at the stationary point, 
is transforming as 



dzidzj d^id^j dzi dzj 
with 



:b.i6) 



dc 



dz 



la — 72;p 1 + zz 



^+« (B.17) 



Consequently if we call Xj the argument of -^ , the similarity property 

BB{x) = e'^BB{z)e-'^ (B.18) 

implies the invariance of the spectrum. The eigenvalue equation can also be written 

A5z + {l-X)D5z = 0, and c.c. (B.19) 

and differs from the naive form related to ( p.8| ) which was written A5z + {D — X)Sz = 0. 
Explicitly, equation (|B.19| ) reads 



A'' 

With the help of ( p.l| ), it can be reexpressed using Szj = Zj^j as 

The 3 zero modes (p = 1) are associated with the 3 eigenvectors 

Szj = iZjd{Arga) or ^j oc i 

5zj = (z| + l)5{Rej) or ^j oc {zj + zj^) (B.22) 

5zj = i{zj — l)6{Im'y) or ^j oc i{zj — z~ ) 
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Let us exhibit one negative eigenvalue in the case of the cuboctahedron. The invariance of 
the system ( p.21|) by the octaedral group reduces the degree of the equations to 3 or 4. In 
the cychc representation of order m = 4, that is taking as a polar axis m the axis having a 
quaternary symmetry, the z/s separate in 4 cosets: y/ii" , ai'^ , a~^i'^ with z^ = 1, 2, 3, 4(mod 
4) belonging to 3 modules 



1, a, a ^, 



a = V2- 1 = tq- 



so that 



12 



f{z) = Hiz - z,) ^ {z^ + l){z' - a'){z' - a-') 



(B.23) 



^12 + 1 _ 33(^8 ^ ^4^ 



We can also choose the ternary representation in 4 cosets 



U 



-U 



bf. b~ r, z^ = 0,l,2mod 3 

X 



b=V3-V2 = tg^ 
with X the angle between ternary and binary axis. This yields the polynomial 

9io = ie + m'+b'){e-b-') 

We switch from ^ to ^ by the rotation 



(B.24) 



(B.25) 



^ zVi + p ^ ii V3- 1 



2 sin 



12 



(B.26) 



where V' is the angle between ternary and quaternary axis. Although we did not want to 
do the complete reduction according to the octahedral group, let us note that this is the 
group of permutations of the 4 diagonals of the cube named 1,2,3,4. The elements of the 
proper group fall in various classes according to the table 



classes 


(1) 


(13) 


(134) 


(13) (24) 


(1234) 


order 


1 


2 


3 


2 


4 


number of elements 


1 


6 


8 


3 


6 


name 


E 




J-i 


/2 


/ 


realization 




z 




z ^ —z 


z -^ iz 
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And the reflection with respect to the center is z ^ —^. 

We are now going to seek the proper modes ( [B.21[ ) of the system corresponding to 
invariant vectors by I (J)- The form ( [B.21D is adapted to Fourier transformation. We hmit 
ourselves to the zero axial moment ^j = ^c where c is the index of the coset of j 



We only need to compute 



E 

jeb 



vi, 



I, I 



Z^Zj 



a, 



a for I m = A 



-6, h-^ for J 



m 



^m^r. 



iZi-ZjY 



m 



(a'" - 6™)2 



I & a 



^B.27) 



where a, b distinguishes here 2 indices of arbitrary cosets. The system ( p.21|) is then 
reduced to 






N-1 



a^-b^Y (|a| + |a-i|)2 



(?a - Pla) = 



(B.28) 



In the case m = A the coefficients are all real because of ( [B.23D -( [B.25P and the system 
decouples in Re^a = or Im^a = according to the sign of p. We are left with one 
equation of second degree in A, after elimination of the root A = in the space of cosets 



v4 = 16 



Z)=ii 



i! 




36x32 

1 


¥ 


36 

-1 


2x242 
31 


36 


2x242 


36x32 


2 


1 


1 







(B.29) 



If we write A' 



9x11 
25 



A, the secular equation 



2-2A' -1 

-1 l-A'-2-5 
-1 2"^ 



-1 

2-5 
1 - A' - 2-5 



(B.30) 



has the three roots A' = 0, 2, 1 — 2 ^ which give the five stable modes 



A 


2-A 


P 


* 

10 
33 

64 
99 


2 

56 
33 

134 
99 


1 

23 
33 

35 
99 
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We would obtain the other modes relative to the moments |/c| = 2, 1 with the ansatz 

e, =ee^'", j = {c,i^), z. = 0,l,2,3(mod4) (B.31) 

but we will now turn to the ternary representation m = 3 limiting ourselves again to the 
total zero moment (states invariant under J). Taking into account ( p.26| ), naming ^1,2,3,4 
the components of the vector ^ relative to the 4 cosets (axial triangles) i, —i, —b, b~^, the 
equations (p.28| ) reads 



-^(6-6) 
7(^-6) 



^(1 + llzy2)(2ei - 6 - ^4) + ^(^1 - P^i) = 
^(1 - llzv^)(26 - 6 - ^4) + ^(6 - P^2) = 



(B.32) 



-4^(6-^4) - ^((26-6-6) - iizv^(6 -6)) + ^(6 -p6) = 

4^(6 - U) - ^((2^4 - 6 - 6) - ll^v^(6 - 6)) + ^(^4 - P^4) = 
Setting 

6 + 6 = ^, 6 + ^4 = ?? 

^(?l-6)=r, 6-^4 = V 

The ternary system (|B.32| ) decouples, as predicted, according to whether the variables 
( |B.33| ) are all real or all imaginary corresponding to 2 values of opposite sign of p and can 
be written 



(B.33) 



11 



35 W ', V c, ; ^ ^ ^2 



■2 35^^ 35 I? 'D 4^32 



1 



iV 



11 



(77 - pf]) 

ie + p'i') 

■ill' - pv') 



(B.34) 



2 34 ' 4 X 33 
which gives assuming that Re{^, ..., rj') = 0, the zero root A = (p = —1), the obvious root 



A = 33^^ and two other solutions of an equation of second degree in /x 



11x3'' 



A 



(p-4)(p-26)-8x 11^ 







(B.35) 



So /U = 15 ± 33. We therefore obtain 7 eigenvalues of which one is negative 



A 


2 - A 


* 




2 


2 
33 


0.0606... 


1.9393... 


26 


= 0.6464... 


1.3535... 


11x32 


2^ 


= -0.2424... 


2.2424... 


11x3 
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We note that the 15 eigenvalues computed (from the 21 non zero ones) are rational numbers. 
Would this be a general property of hessian matrices of equilibrium configuration on the 
sphere? We should first check this on the icosahedron. 

It is interesting to know the direction of the unstable ternary mode A = ^ = 
-0.2424... and also that of the stable mode A = ^ = 0.0606.... With the help of ( WM ) 
we obtain for the unstable ternary mode the proportion 

Im{^ : ^' : 1] : 1]') = 1 : -V2 : -3 : (B.36) 



and for the less stable only Imr]' is non zero. This gives according to (p.33|) for the unstable 
mode 

^3 = ^4 = -1.5Z 

(B.37) 
^1 = O.bi - 0.707... 6 = 0.5z + 0.707... 

and for the stable mode 

6=6 = 0, ^3 = -^4 = ^ (B.38) 

The elementary displacements dz in C are such that Sz = ^z according to ( [B.21| ). To 
interpret the result ( p.35| ) we see that it is possible, by a global rotation around the 
ternary axis {Szj/zj = 1.5z), to keep fixed the two cosets —b and b~^ made of the two 
opposite triangles of the cuboctahedron centered on the ternary axis. Finally this gives us 

the displacements 

dzs = dz4 = 

1 

^72 (B.39) 



which describe a deformation of the plane hexagon of the cuboctahedron made of i) a 
rotation of the hexagon (— = 2i) ii) a deformation of the latter by dedoubling of the plane 
of the two cosets i and — i The value of the angle {Arctg{l/{2y/2)) and the direction 
indicated on the figure suggests that this unstable ternary mode describes well the begin- 
ning of the displacement of the cuboctahedron towards the stable icosahedron, the two 
opposite triangles staying the same to first order. The less stable mode ( [B.38[ ) is related to 
the vibration of torsion (opposite rotations) of only these two triangles, the plane hexagon 
remaining fixed. 

Fig. 6: Going from the cuboctahedron to the icosahedron 
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6zs = 


6Z4 = 


-0 


6zi/z 


L=2Z 


A 


SZ2 

Z2 


--2i + 


1 

V2 



Appendix C. Shiota's Formula 



It is possible to write a formula for A in terms of the traces ti = X]i<i<Ar ^l- Noting 



that 



2s 



n ( 

l<i<j<N 



Zi - Z 



\2s 



I ^Jy 



is (-iyN{N-i)/2 ^ij^gg ^Yie coefficient of e^(^-i) in 



- e(Zi - Zj^ 



l<i<N 
l<j<N 



where e is a formal scalar parameter close to zero. Let us compute this quantity 

\ 



n = exp Y^ log(l - e( 



l<i<N 
l<j<N 



exp 



exp 



exp 



n=l l<i<iV 

l<j<iV 



Zi Zj ^ 



Zi Zj J 



\ 

s 



°° ^n ^--^ /^„\ 



Tl ** ^ \ / 

n=l 1=0 ^ ^ l<i<N 

i<j<N 



I 



oo „ ns 



n '■ — ' V / 



E^E T(-i)'*"-'*' 



Taylor expanding the exponential we have 



(C.l) 



(C.2) 



(C.3) 



n-EM-E^E Tm)V,-,. 

A;=0 \ n=l 1=0 ^ ■' > 

Taking (— l)''^(^~i)/2 times the coefficient of e^(^-i) we obtain 



Z^ 11 9^ ^ 



^2s ^ (^_i)sJV(7V-l)/2y^ 



A;! 



n,s 



2?i,- -^^ V / 

k>\ «i>i,...,»fe>i i=l 1=0 

ni + ...+nj,=JV{JV-l) 



(C.4) 



(-l)^t^^,_z tz (C.5) 



(summing over k < N{N — 1) suffices since the ni > 1 must add up to A^(A^ — 1) ), or 



X 



E 



/_l\ki+k2 + 



fcl ,1^2, ■■■>0 
ki+2k2+3k^ + ... = N(N-l) 



ki\k^\ 



n 



1 ^^^ ns 



, 71 -^-^ V / 

n>l \ 1=0 



( — 1) t„s_i t; 



(C.6) 
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Appendix D. Decomposition in terms of characters for small A^ 

The decomposition of the square of the Vandermonde determinant in terms of char- 
acters can be written in the following form, 



A^ = J](i - T,)\l - nn+if ■ ■ ■ (1 - n...TN-i] 



c/l0,3,...,3(Af-2),3(Af-l) (D.l) 

adm 



where the operators Ti are acting on the labels of the characters at position i shifting li-i 
upward by one and li downward by one. Doing the decomposition explicitly we find 

N=2 (2 terms) 

A^ ^-3chi,2 + cho,3. (D.2) 



N=3 (5 terms) 



A^ = -15c/l2,3,4 + 6c/li,3,5 - 3c/lo,4,5 

- 3c/li 2,6 + cho- 



[D.3) 



N=4 (16 terms) 



A = 105c/l3,4,5,6 - 45c/l2,4,5,7 " 6c/l2,3,6,7 + 27c/li,4,6,7 
- 15c/lo,5,6,7 + 27c/l2,3,5,8 " 9c/li,4,5,8 " 12c/li,3,6,8 
+ 6c/lo,4,6,8 + 9c/li,2,7,8 " 3c/lo,3,7,8 " 15c/l2, 3,4,9 
+ 6c/li,3,5,9 — 3c/lo,4,5,9 — 3c/li,2,6,9 + c/lQ, 3,6,9- 
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(D.4) 



(D.5) 



N=5 (59 terms) 

A^ = 945c/l4,5,6,7,8 - 420c/l3,5,6,7,9 - 75c/l3,4,6,8,9 
+ 270c/l2,5,6,8,9 + 45c/l2,4,7,8,9 - 180c/li,5,7,8,9 
+ 105c/lo,6,7,8,9 + 270c/l3,4,6,7,10 - 90c/l2,5,6,7,10 
+ 45c/l3,4,5,8,10 - 144c/l2,4,6,8,10 + 72c/li,5,6,8,10 

— 18c/l2,3,7,8,10 + 81c/li,4,7,8,10 " 45c/lo,5,7,8,10 

— 18c/l2,4,5,9,10 + lllc/i2,3,6,9,10 — 27c/li,4,6,9,10 

— 6c/lo,5,6,9,10 ~ 54c/li^3j^9^io + 27c/lo,4,7,9,10 
+ 45c/li,2,8,9,10 - 15c/lo,3,8,9,10 " 180c/l3,4,5,7,ll 
+ 72c/l2,4,6,7,ll - 36c/li,5,6,7,ll + 81c/l2,4,5,8,ll 

— 27c/l2,3,6,8,ll " 36c/li,4,6,8,ll + 27c/lo,5,6,8,ll 
+ 18c/li,3,7,8,ll - 9c/lo,4,7,8,ll " 54c/l2,3,5,9,ll 
+ 18c/li^4^5^9^ii + 24c/li^3^6^9^ii — 12c/lo,4,6,9,ll 

— 18c/li,2,7,9,ll + 6c/lo,3,7,9,ll + 45c/l2, 3,4, 10,11 

— 18c/li^3^5^io,ii + 9c/lo,4,5,10,ll + 9c/li^2,6,10,ll 

— 3c/lo,3,6,10,ll + 105c/l3,4,5,6,12 " 45c/l2,4,5,7,12 

— 6c/l2,3,6,7,12 + 27c/li,4,6,7,12 - 15c/lo,5,6,7,12 
+ 27c/l2,3,5,8,12 - 9c/li,4,5,8,12 " 12c/li,3,6,8,12 
+ 6c/lo,4,6,8,12 + 9c/li,2,7,8,12 — 3c/lo,3,7,8,12 

— 15c/l2,3,4,9,12 + 6c/li,3,5,9,i2 — 3c/lo,4,5,9,12 

— 3c/li^2,6,9,12 + c/lo,3,6,9,12- 

If we label the coefficients by C{Ti...Tj) the following simple rules can be found directly 
by looking at ( p.l|) , 

C(riri+i...ri+„) = C(rjrj+i...rj+„) (D.6) 

C{r,...r,^^) = -2C(r,+i...r,+p) (D.7) 

We also read from these examples the following list of coefficients, 

Cin) = -3 (D.8) 

C{rf) = 3 (D.9) 
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C{n...n^^) = -?>{-2r (D.io) 



C(r/...r4J=3(4--9n) (D.ll) 



27 99 1 



C(r.rii...r^Vi) = -3(2^-^ " -^ + - + 7(-5)"-') (D-12) 



T 



C(r,rf+i...r2+^_2r,+^_i) = 162 + (24 - -n)(-5)"-4 + 3 22-4 _ __^ (j3_^3) 



For higher powers of the Vandermonde determinant we found the foUowing decompo- 
sitions in terms of characters. 

Fourth power 

N=2 (3 terms) 

A^ = c/io,5 - 5c/ii,4 + 10c/i2,3 (D-14) 



N=3 (13 terms) 



A"^ = c/lo,5,10 - 5c/li,4,10 - 5c/lo,6,9 + 10c/l2,3,10 + 10c/lo,7,8 " 15c/l2,5,8 
+ 20c/li,5,9 - 25c/li,6,8 - 25c/l2,4,9 + 100c/l3,4,8 + 100c/l2,6,7 (D-15) 

- 160c/l3,5,7 + 280c/l4,5,6 
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N=4 (76 terms) 

A = c/lo,5,10,15 ~ 5c/lo,5,ll,14 + 10c/lo,5,12,13 

- 5c/lo,6,9,15 + 20c/lo,6,10,14 " 25c/lo,6,ll,13 + 10c/lo,7,8,15 

- 25c/lo,7,9,14 - 15c/lo,7,10,13 + 100c/lo,7,ll,12 + 100c/lo,8,9,13 
+ 280c/lo,9,10,ll ~ 160c/lo,8,10,12 - 5c/li,4,10,15 + 25c/li,4,ll,14 

- 50c/li,4,12,13 + 20c/li,5,9,i5 - 80c/li,5,io,l4 + 100c/li,5,ii,i3 

- 25c/li,6,8,15 + 50c/li,6,9,14 + 100c/li,6, 10,13 - 375c/li,6,ll,12 
+ 75c/li,7,8,14 - 300c/li,7,9,i3 + 375c/li,7,io,i2 + 300c/li,8,9,12 

- 600c/li,8,10,ll + 10c/l2,3, 10,15 - 50c/l2,3,ll,14 + 100c/l2,3,12,13 

- 25c/l2,4,9,15 + 100c/l2,4,10,14 - 125c/l2,4,ll,13 " 15c/l2,5,8,15 

+ 100c/l2,5,9,14 - 290c/l2,5,10,13 + 475c/l2,5,ll,12 + 100c/l2,6,7,15 

- 300c/l2,6,8,14 + 200c/l2,6,9,13 + 100c/l2,6,10,12 + 825c/l2,7,8,13 

- 1125c/l2,7,9,12 + 150c/l2, 7,10,11 + 1200c/l2,8,9,ll + 100c/l3,4,8,15 

- 375c/l3,4,9,14 + 475c/l3,4,10,13 - 250c/l3,4,ll,12 - 160c/l3,5,7,i5 
+ 375c/l3,5,8,14 + 100c/l3,5,9,i3 - 685c/l3,5,io,l2 + 300c/l3,6,7,14 

- 1125c/l3,6,8,13 + 1425c/l3,6,9,12 " 750c/l3,6, 10,11 - 200c/l3,7,8,12 
+ 800c/l3,7,9,ll - 4400c/l3,8,9,10 + 280c/l4,5,6,15 - 600c/l4,5,7,14 
+ 150c/l4,5,8,13 - 750c/l4,5,9,i2 + 3180c/l4,5,io,ii + 1200c/l4,6,7,13 
+ 800c/l4,6,8,12 - 3800c/l4,6,9,ll - 1200c/l4,7,8,ll + 6600c/l4,7,9,io 

- 4400c/l5,6,7,12 + 6600c/l5,6,8,ll - 880c/l5,6,9,10 - 9240c/l5,7,8,10 
+ 15400c/l6,7,8,9 



Sixth power 

N=2 (4 terms) 



(D.16) 



A^ = c/lo,7 - 7c/li,6 + 21c/l2,5 - 35c/i3,4 (D-17) 
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N=3 (25 terms) 

A*" = c/lo,7,14 - 7c/lo,8,13 + 21c/lo,9,12 
- 35c/lo,10,ll - 7c/li,6,14 + 42c/li,7,i3 - 98c/li,8,12 
+ 98c/li,9,ii + 21c/l2,5,14 - 98c/l2,6,13 + 119c/l2,7,12 

+ 147c/i2,8,ii - 539c/i2,9,io - 35c/i3,4,i4 + 98c/i3,5,i3 (D-18) 

+ 147c/l3,6,12 - 868c/l3,7,ii + 1078c/l3,8,10 - 539c/l4,5,12 

+ 1078c/i4^6^ii + 231c/i4j^io — 2695c/i4^8^9 ~ 2695c/i5^6,io 
+ 3850c/i5 7 9 — 5775c/i6 7 8 



Eighth power 

N==2 (5 terms) 

A^ = c/io,9 — Qchi^s + 36c/i2,7 

- 84c/l3,6 + 126c/l4,5 

N=3 (41 terms) 



(D.19) 



A^ = c/lo,9,18 - 9c/lo,10,17 + 36c/lo,ll,16 

- 84c/lo,12,15 + 126c/lo,13,14 - 9c/li,8,18 + 72c/li,9,i7 

- 243c/ii,io,i6 + 432c/ii,ii,i5 - 378c/ii, 12,14 + 36c/i2,7,i8 

- 243c/l2,8,17 + 603c/l2,9,16 - 432c/l2,10,15 - 972c/l2,ll,14 
+ 2646c/l2,12,13 - 84c/l3,6,18 + 432c/l3,7,17 - 432c/l3,8,16 

- 1776c/i3,9,i5 + 5724c/i3,io,i4 - 6048c/i3,ii,i3 + 126c/i4,5,i8 (D-20) 

- 378c/l4,6,17 — 972c/l4,7,i6 + 5724c/l4,8,15 — 7650c/l4,9,i4 

- 3402c/i4,io,i3 + 19656c/i4,ii,i2 + 2646c/i5,6,i6 ~ 6048c/i5,7,i5 

- 3402c/i5,8,i4 + 28224c/i5,9,i3 - 34398c/i5,io,i2 + 19656c/i6,7,i4 

- 34398c/i6,8,i3 - 3822c/i6,9,i2 + 68796c/i6,io,ii + 68796c/i7,8,i2 

- 91728c/i7 9 11 + 126126c/i8 9 10 



CoefRcients C/\ for the vertices of the admissible polytope. 

The admissible monomials in the expansion (D.l) are of the form t^^ ...Tj^Zi , where 
the (no = 0, ni, ..., un-i, un = 0) are subject to 

2s + Ui+i + Ui-i - 2ni > 0; m > 0; i = 1, ..., A^ - 1. (D-21) 
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The equations ( p.21| ) define a polytope or IR ~ , which is a certain deformation of an 
hypercube. It can be shown [|12[ that its 2-^"^ vertices are in one to one correspondence 
with strictly increasing sequences of integers: 

zo = < zi < 22... <ik < ik+1 = N 

denoted [0,ii, ...,ik,N] here. The corresponding points read 

Ui^ s{i-ik){ik+i-i), Vi G [zfc,Zfc+i], / = 1, ..., /c - 1. (D.22) 

These points are easily identified as the only ones which vanish at i = ii, ...,ik (saturate 
the second inequality of ( p.21|) at these points), and have Ui^i — 7ii_i — 2?!^ + 2s = 
(saturate the first inequality of (|D.21|) ) between those. 

The farthest vertex from the origin corresponds to [0,A^], and rii = si{N — i) for 
z = 1, ..., A^— 1. It yields the most compact character chs(N-i),s{N-i)+i,...,{s+i){N-i) with 
coefficient already given in (5.16) 

Ms) _ [(g + l)A^]! , .yN(N-l)/2 (j^2'^) 

Thanks to the factorization property 5 (eqn ( [4.43| ) ) the coefficient pertaining to the vertex 

[0,ii, ...,ik,N] reads 

k 






1=0 



Appendix E. More on the number of terms in the expansion 
Computing the volume 6]y = a}^ = Vn+i 



If we reorganize the inequalities ( |8.1|) for s = 1 , we find that 

nN-2 



sup{0, 2nN-2 - riN-s - 2) < un-i < 1 + 
SUp{0, 2nN-S - "IT-N-A - 2) < n7v_2 < 2 + 



2 

2nN-3 



sup{0, 2n2 - ni - 2) < na < iV - 3 + ^^ _^^"' 

sup{0, 2ni-2) <n2<N -2 + ^— — ^^ 
< ni < A^-1 
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Let us define (j)N-i-k as (/)n-i = 1, and 



{k + l)nN-k-2 



k+- 



fe + i 



(pN-l-k = 0{k+l^ — — HN-l-k) / dnN-k(pN-k, 

f^ + ^ JsMp(0,2njv-i-fc-njv-fc-2-2) 

where 

e{x) = 1 if X > 

— otlierwise. 
For instance, 

(/)7V-2 = 0{2 H nN-2)inf{l H ^, -(2 H nAr_2)). 

Tlie volume of the poly tope 7r]y is 

Vat = (pQ. 

Let us rewrite 

(t)N-2 = 2(2H ^ nN-2)0{2-{ nAr_2)-2(H nN-2)0{U nN-2): 



and introduce the notation 



0) . , JUk-i 

""k =3 + -qry - ^k 



(3^1 = J ~ f^i and Xq = j) then 



?Ar_2 — r,''^N-2^\''^N-2) ^.X j^_2U\X j^_2) 



Suppose we know that (t)N-k has the form 



^N-k=J2pf~'^^^^N-k)0{4lu) (E.l) 



where the p^^ ^-*'s are some polynomials of degree k — 1. Then by definition, 

, . fe"JV-l-fc 

<^Ar-fc-l = 0{x)^tk_i) I duN-k (pN-k 

-'sMp(0,2njv-i-fc-njv-fe-2-2) 

kn ]• 



j = l "'s^ip(0,2njv_i_fe-njv_fc_2-2) 



.■_i J suv(0.2nAT--i -h — nM-k.-o — 2) 
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ii) 



which becomes after the changes of variables riN-k -^ ^N-k 

k 
-n(Ak+^) \\^Q(,^f(A I J^N-l-k j + 2 (j+i) 



)N- 



J = l 






j + l 'jH 



du p- (u) 



Note that 






X 



AT-fc- 



-1-0' + 



JUN-l-k 
J + 1 



2x 



(1) 

Af-fc-l' 



SO that 



^N-k-l 



E[i - o(4lk-im4-ti) f'^'^"-'-' du pf-'\u) 






VI I "iV-l-fe \ 



du p\ ' {u) 



^ "Jil+^^^^J^) 






±)+2x(l) 



du p- (u) 



iV-fe-l 



-E''( 



X 



(i+i) ^ 

AT-fc-l. 



j + 2^(j + l) 

^■ + ^^-'=-1 (fc-1), ^ 



where we used the inequahty x]y_^_j^ < x]^_^_j^ for j = 1, ..., /c to rewrite 



The form of (pN-k-i will be that of the recursion hypothesis ( JE.ID iff 






(E.2) 



is independent on y, and a polynomial of x only. Then we have, for j = 2, ...,k + 1 

■ i+ix 

d.ii. 1^ 




p) '[x) 



du p)_^ {u)- 



Thanks to these relations, we can write p\ as 

fc+i 



Sk), 



Pi \x) = Y1 

j=2 






i + i 



i j + i' 
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Let us prove by recursion that this expression is independent on y. Suppose the 
property true for /c — 1, then in the above we can use it to rewrite ( [E.2|) 



p['\x) = j2l 



du p- (u) 



=i^3a+jh)+2x 



by replacing 



k r 



P 



J=2 

vahd for any z. We get 

J=2 



J + l J ■' J + l 3 J + 1 



i(l + 7^)+2x 



du p\ (u) 



+ 



£L^'""<f:^<-f"-^'""'T^'-^^ 



J + 1 J ■" J + 1 J J + 

We take the range of the integrations to be [2x, 0] by suitable changes (shifts) of variables, 
to get 

r-O 1^ 



pf\x) 






,(fc-l) 



y 



p'; ^>{u+j{i + -^)+ 



+pr"(^(i + -)) -pr"(^(i + ^) + ^(u + 1 + f )) 

■^ J + l J -^ j + l jj + l 2 

But this is valid for any z, let us take z = —u, then the first and last term in the sum 
cancel out exactly, leaving us with 

k „o 



p[^\x) 



t^ i2x ^ J + l J 



clearly independent of y. This completes the proof of the general formula 

fc+i 






with the recursions 

p[ (x) = 1 



<M, 



3 + 1. 



p';\x)= dup)_,>iu) i = 2,3,...,A; + l 

Jo 

fc+i 



.(fc)/ 



Pi \x) = Y1 

j=2 



j ^ i + 1 



i + i 



i j + i' 
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where the last expression is independent of y. 



The volume reads 



TV 






(fc)/ 



The expression for p\ [x) gives for y = 



fc+i 

Pi \x) =Y^ 
i=2 



.(^), 






hence for x — 1, 



k+l k+l ., . ^x 



The proof of independence of y gives us another expression for jy^ 



(k) 



1 Z 7 ~r J- 7 J- 



J=3 



J 



J + 



F^)] 



The first few p's are listed below. 
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A; = pf^ = 1 
k = l p^^^ = -X 

p\ = —2x 
k = 2 pf^ = ^x^ 

(2) 9 2 
P2 = -^^ 

p\ = x"^ — 2x 
^ o (3) 125 3 

(3) 16 3 
^3 g 

(3) 9 3 9 2 
p\ = -X -X 

pf^ = - x^ +2x'^ -4x 
9 



9 

— X 
20 



,4 



k = 4 pI^^ = 

"576 



4 625 4 

p\ = X 



(4) 84 lU 3 

P3 = -x — :^a; 



16 
9"" Y 

f4) 949392 
P9 = X + -X X 



32 4 2" 

(4) 1 4 232 32 

p) = — X X + ox X 

^^ 36 3 3 

^^ 86400 

(5) 27 5 
^^ 50 

(5) 625 5 625 4 

p^ = X X 

^^ 1152 576 

(5) 32 c 16 4 32 3 

p\ = X H X X 

^^ 135 9 9 

(5) 27 c 27 4 45 q ,„ 9 
p^^^ = -—x^ - —x^ + —x^ - 12x2 



640 32" 
^ 5 1 _ 

450 9 3 3 
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(5) 1 5 1 4 o 3 44 2 100 

p) = X -\ — X — 2x -\ X X 



which lead to the volumes 



N 


Volume Vn 


2 


1 


3 


2 


4 


16 
3 




50 


5 


3 


6 


288 


5 


7 


9604 
45 


8 


262144 



315 



It is easy to read the actual exact result 



(AT) _ .(N) _ . . _r,N 



{N + 1 



,N-2 



N\ 



(E.3) 



hence the dilated poly tope n^;^^^ has volume 



'(s) _AN)N 



^ATi = <^«" = (2«) 



N 



{N+1 



,N-2 



m 



(E.4) 



We guess formula (E.3) on the basis of the previous calculations. 
Computing 6)v-i ~ ^n-i 



Let us now concentrate on the subleading coefficient ajy^ = b}^_i- By definition, we 



have 



AN) N-i_^l(_d_ d^.yis) ., 



k=l 



2 dhk dek 



(E.5) 



hi=ei=Q 



Let us write the deformed volume 

N 



Vl^%^{h,-e^)= f d^u\\e{uk+ek)9{-2uk + Uk+i+Uk-i+2s + hk). (E.6) 

^^(x) = (5(x). Therefore we have two contributions to ([E.5[) , the one coming from e 
derivatives 



r{s) 



de^V^li{K]ei] 



hi=ei- 



N 

= / d^uS{uk)Y\d{uj)Y\d{-2ui + ui+i+ui_i + 2s) 

Jueu'^ /4 tl 



^k ^N-k 



85 



and the one coming from h derivatives 






k=i 



n+i 



is) 



Putting all contributions together, we find 



AT 1 ^ 

7 AT-l AT-l ^^ ir(s) , -*- V^ t/(s)t/('5) 



4s ''+' 2 

A;=0 



N-k- 



The latter sum turns out to have a simple form 



^(s)^As) _ (2g) 



E ^^''^N-^ = (^3T)! (^ + l)''"'^^ + ^)' 



which leads to 



nAT-l 

^''-^i = ^''-^i = (]^3T)T^(^ + 1)'^"'(6 + (iV + l)(iV + 2)). (E.7) 

Applying this to A^ = 2, 3, 4, 5, 6, 7, 8, 9 yields respectively 2, f, 24, 96, ^, ^^, ^^|p, 
2320000 gQ ^g gg^ ^Yl the subleading coefficients in the list ( p.2|) . Note that the general 
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expression fails to reproduce the corresponding A^ = 1 term, due to the different definition 
of the deformation for A^ = 1. We also get the subleading coefficients in the expression of 
the A's in terms of the a's ( |8.3| ), namely the coefficient of ctn-i in ^n+i rcads 

2^-2(Ar + l)^-4[6 + A^(A^ + 1)(2A^ + 1)], (E.8) 

which yields for A^ = 2, 3, 4, 5, 6, 7, 8, 9, 10 respectively 4, 45, 744, 16128, 432768, 13860864, 
516481920, 21964800000, 1050351430656. 

Computing bj^_2 and a]Y_2 

The next to leading coefficient a^_2 ^^ trickier to obtain, because it differs from bN-2 
by a highly non-trivial quantity. The latter involves the two-dimensional cone invariant 
T2 computed in [10] and function of the Dedekind sum of the two integers defining the 
(rational) cone. However in our very particular case, we are able to guess a simple answer. 
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The first task is to compute bN-2 then to proceed and compare its value to that of aN-2 
which we read from ( p.2|) . Thanks to the definition (|8.8| ), we find a closed formula for 67V-2 



bN-2s''-' 



i i<j 



+ zE^^^^^ 



«>j 



|/^_|_l(,'ii5 £i 



hi=ei=0 



We use the fact that 



Y.^r..V 



d 



d{2s 



-V 






to rewrite 



bN-2S 



N-2 



+ 



12 

8ds 



^E(9i-H) + iE*.*.+ 



E* 



r(s) 



V^ ' (h- e- 



i<j 



hi=ei=0 



J_d^ 

?,2ds' 



Each of the above terms is easily derived by using the integral expression for the deformed 
volume ( [lii.6| ). We find respectively 

^d^^d,.VN+l= J] VpVqVN-p-q 



i<j 



d 



p,q>l 
p + q<N-l 



ds- 



-Vi 



N+1 



N{N-l)VN+i/i^ 

N-l 



-rJ2d,,VN+i 



N-l 



ds 



J2 ^p^N- 



p=i 



The only difficult part involves the double derivatives. For e derivatives, we have 

dlVN+i = f d^ud'{u,)Y\9{uj)Y\9{-2uj + Uj+i + w,_i + 2s) 
JueTR'^ V; , 



Rf^ 



du5{ui) '^6{uj)(pi{uj), 



j^i 
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where 

^i{uj) = 5{ui + Ui+2 + 2s - 2ui+i) JJ^ 9{-2uj + Wj+i + w-,_i + 2s) 

+ 5(wi + Ui-2 + 2s — 2ui-i) TT 6{—2uj + Wj+i + tij-i + 2s) 
- 2(5(ui+i + Wj-i + 2s - 2ui) Yld{-2uj + Uj+i + Uj-i + 2s). 

The last term does not contribute because it multiphes d{ui)9{ui-i)9{uij^i) and for s > 
the various constraints are incompatible. So we are finally left with 

N-l 
p=l 

Using again the integral definition (|E.6| ), we find 

1 '~^ 

so that 

p + q<N-l 

For h derivatives, we finda 

^ (2s)^-^(iV + 1)^-3 

i 

Finally, we use 

.(.).....) 3 (2s)^-^(iV + l) ^-^ 



V y(^V(^Vr^ = -^^ \fUL±I \N^ + 15N + 741 

Z^ ^P ^q ^N-p-q 4 (N-2)\ -r±oiv -r(^j, 

p,q>l 
-q<N-l 

SO that bN-2 reads 



(A^-2) 

p,q>l ^ ' 

p + q<N-l 



>iv-2 = ^,\^_^ , [3iV^ + 17Ar3 + 72N^ + lUN + 266]. (E.9) 



Now comparing this to the values of aN-2 that we read off from (p?2|), we conjecture 
the following relation 

a^.^s""-' = 6Ar-2«^-' + ^4'2i, (E.IO) 



See next section for a detailed proof 



leading to 

2Af-2 

a^_2 = [{N+l)^-^{3N^ + 17N^ + 72N^ + U4N+266)+6{N-l)^-^]. (E.ll) 

We check that for A^ = 4,5,6,7,8,9,10, we get aN-2 = f , ^, ^, ^^, ^^g^, 
39666608 ^ 205427098 p^^^g ^y^^^ ^j^g g^^^^^g formula fails to reproduce the A^ = 2, 3 cases. 

This, we suspect, is related to the difference in the definition of the deformation for the 
boundary case A^ = 1. 

This leads to the coefficient of o"Ar_2 in A^^j^^ of (|8.3| ) 

2iV-2 

— j— [(A^+1)^-'^[12A^*^+16A^^-25A^^+21A^3+16A^2-100A^+114]+6(A^-1)^-^] . (E.12) 

For A^ = 4, 5, 6, 7, 8, 9, 10, we get respectively 404, 12481, 437776, 17367872, 773038912,38261688576, 
2088303502080. 

Computing the polynomials B^^s) 

The computation of the polynomials i?Ar_i(s) resembles very much that of the volume 
Vn{s) of previous appendix. The only difference is that the defining relations for the 
polytope are decorated with the parameters hi and e^, i = 1,2,...,A^— 1. For simplicity, 
let us first compute the deformed volume WN{hi, e^) defined by the relations 

-2efc < 2nk < Uk-i + Uk+i + hu- 

The desired volume is 

V^^'\h,,ti) = WN{K + 2s,ti). 

Let us again introduce (j)N-i{hi, e^) = 1 and 07v-fc-i(^i, Ci) by the recursion 

fc+i 
(t)N-i-k{hi,ei;x = UN-i-k) = 0(^jhN-j + {k + l)nN-k - {k + 2)x)x 

du (j)N-k{hi,ei;u). 

'sup( — ejv_fe,2njv_i_fc— njv_2-fe — /ijv-i-fe) 

With these definitions, the volume W reads 

WN{hi,ei) = (f)o{hi,ei]0). 
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Just as in the non-deformed case, the (/)'s turn out to have a simple polynomial form 



(/i„ q; x) = Y,Pj '^ (4-fc(^- e^; x))0{x%l,{h,, e,; x)), (E.13) 



where the p- (u), j = 1, 2, ..,1 + 1 are some polynomials of u, with coefficients themselves 
polynomial in hi and e^, and we used some reduced variables 



1 

XN_,^{hi,ei]x) = — — (^^ rhN-k+j-r - ^N-k+j) + -—-nN-k-i -x 

^ '^ r=l ^ "*" 



with the convention that cn = 0. The proof of ( E.13|) is identical to that of the undeformed 



case and leads to the recursion relations between the p's 



pf\x)=l dupf_\'\u) j = 2,3,...,k + l, (E.14) 

Jo 



and 



Wr \ sr^ ik)fYli=irhN-r+jeN-k-eN-k+j-i + (j-l)y 

p\ (u) = 2^p)' y-^ j^^ 



p 



J 



-j-i 



(E.15) 



(fc) fT.r=l rhN-r+j(^N-k - eN-k+j-1 + {j - l)y + 2j 



u 



j + l 



independent on y. 

We list the first few deformed p's below 
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k = pl {u 

k = 1 P2 (u 

Pl [u 

k = 2 pf\u 

(2)/ 

p\ (u 
k = 3 P4 (w 

(3)/ 
(3)/ 

P2 (« 



1 

3 



-w 



-2m 



3" 



-r 



V? - [/lAT-l + 2eAr_i - eAr_2]w 

125 o 

144" 

16 o 

u-^ 

9 

9 . 



.[/. 



N-\ 



+ 2eAr_i -eAr_2]w 



--W + -[/lAr_2 +2/lAr_i -f 3eAr_i -eAT-sJM 

- -\^N-\ +h%-2 +4:hN-lhN-2 + hN-2{GeN-l - 2eAr_3) 
D 

+ hN_i{6eN-2 - 4eAr_3) + e%_^ - 3(e|r_i + e%_2) + 12eN-ieN-2 - QeN-ieN-3]u 

We recover the undeformed p's (at s = 1) by taking /i^ = 2, V z. A few remarks are in 
order in view of this hst. 
(i) The coefficients of the polynomials p- have a simple structure for j > 1, due to the 



recursive definition ( [bil.l4| ). Let us denote by ctm' the coefficients of pi 



(i,fc) 



Ak) 



pf\^) = Y.^'^ 



k) m 



771 = 1 



(k), 



where we used the definition ( [E.15| ), to include the fact that pi (0) = 0. Then the 
coefficients am of p^ satisfy the recursion 



a 



U,k) 
m+1 



1 _^l±l-)m+l^U-l,k-l) ^^1^2 /c-1, 



m + 1 j 



which implies 

,(k+l,k) _ 
n 



(fc + 2)^ (1,0) 
^^■\k + l)\k\^'' 



6{rn — j)0{k — m 



C7±ir>i-j + i)! fi,fc-,+i) 



j!m!2"^~-?'+i 
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a 



m— j + l 



i = 2,3,...,/c. 



So we have only to compute the coefficients aL' oi pl . For instance, due to ctg ' 

1, a]^ = —2, 02 = 1, «! = — (/lAT-i + 2eAr_i — eN-2), we have 



j = k — 1 : a 



{k + iy.ki 

{k,k) _ {k + lf 

ik\y 
ik-i,k) _ {ky 



j - ^ ■ ^k - (u\\2 



k 



k\{k -l)\2\ 



(k-l,k) {kf ^ /r , o \ 

«fc-i = ~ 2\((k - 1)!)2 ^^~^ ^^"^ ~ ^^-'^'- 

The coefficients of p\ satisfy the following recursion relation inherited from the def- 
inition (|El^ 



m'.ik — 'm)\{k + 1)\ ^-^ 



T^{k^\) 
m\(k — m)!(/c + _,. 

(E.16) 
for m = 1, 2, ..., /c. As a consequence, we derive the leading coefficient a\' of p^ 



(E.17) 



°'^' ' A;!(A; + 1)! ^ -^ k\j\ "'^-J+i 



-2 



ifc 



(A;!)2 ' 
thanks to the identity 

E(x-jr(-iy(^^^)=o, (E.18) 



\k 



J=0 

valid for any x, taken at a; = 

k 

j=i \ J / 

(ii) The dependence on hi and e^ of the coefficients am is very peculiar too. Actually, 
the coefficients am ar( 
a consequence of ( [E.16|) . 



the coefficients am are functions of /iat-i, ..., hN-k+i and ctv-i, ..., eAf-fc only, as 
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Finally the volume reads 



TV 






'' ' i+1 



This provides us with a powerful recursive scheme for computing the deformed volumes. 

As an example let us derive exactly the coefficient of/i^ ~^, i = 1,2,...,A^— 1, in Wn. 
The above remark (i) enables to write the leading coefficient of the p's as 

(j,fc) ^ (j + l)^(A:-j + l)! (i,fc-,+i) 

fc j|^l2A;-i+l "fc-J + 1 

^ (_l)fc-i+i_ (-^' + 1) 



j\k\{k-j + l)\ 

(k) 

The contribution to Wn of any subleading piece, say of degree m < k, of p- involves, 
thanks to remark (ii) a coefficient which is only a function of /iAr_i,...,/iAr-fc+j-i, and which 

multiplies a power 

i-i 

r=l 

but A^ — /c+j — 1 >j — 1, therefore such a term does not contribute to any /i^ ~ term of 
Wn- The only contributions come from the leading pieces of the p's 



N A3,N-l) j-l 
7 = 1 ^•' ' r=\ 



and we get the coefficient of h- in Wn 



N N 



Eo-»)-'4^-"-E '^'"^"'''"""^ 



j=i+l i=i+ 

JV-i-l 



^.f^^ j!(A-l)!(A-j)! 



N\(N _,. 



j^.gi-m^-^r-f^) 



N\(N _,. 
^ ' j=o 
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The last identity expresses the symmetry of Wn under the interchange hi <-^ /iat-^ and is 
a consequence of the identity (|E.18|) , taken at x = i and k = N — 1 



N / s 



i=o 
FinaUy we have at leading order in the /i^'s 

N-l ,N-1 i-l /^\ 

This leads to the following contribution of the (A^ — l)th derivatives of Wn w.r.t. the /I's 

N-l N-li-1 . . 

1 

~ N' 
Using this recursive method, we computed the first few -Byv's 

Bois) = 1 

Bi{s) = s + l 

19 

B2{s) = 2s^ + 2s+- 
6b 

^ , , 16 3 13 2 8 107 

Bs(s) = — s"^ H s^ + -s-\ 

^^ ^ 3 2 3 288 

„ , , 50 4 o 77 2 109 641 



3 6 36 2400 

r. r ^ 288 5 ^^4 377 o 1453 2 509 51103 

^ ^ 5 6 72 160 259200 

„ , , 9604 R 6076 . 5641 4 6821 o 20111 o 8663 1897879 



45 15 18 54 720 2700 12700800 
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